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-J  Floodplain  management  has  become  a  topic  of  high  national  priority  in  re¬ 
cent  years.  Computer  simulation  program?,  such  as  HEC-1  and  KBC-T ,  have  great 
increased  the  efficiency  with  which  calculations  can  be  performed  once  input 
data  are  properly  formatted.  The  objective  of  this  study  was  to  develop  systei 
atic  procedures  for  the  conduct  of  river  basin  floodplain  studies.  Guidelines 
are  presented  for  standardising  watershed  subdivision  and  the  collecting,  pro¬ 
cessing,  and  storing  of  stream  cross  section,  elevation,  soils,  land  use,  and 
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other  information  in  computer  data  bases.  Methodologies  for  accessing  data 
bases  and  generating  input  decks  for  computer  programs  HEC-1  and  HEC-2  are 
presented.  Also  discussed  are  procedures  for  interfacing  these  two  programs 
such  that  manual  Input  requirements  are  minimized.  Other  topics  addressed 
include  the  coding  of  bridges,  the  effect  of  reach  length  on  computed  water 
surface  elevation,  the  areal  and  temporal  distributions  of  rainfall,  model 
calibration,  the  autoaiated  mapping  of  floods,  and  flood  damage  estimation. 
Three  software  routines  and  some  pertinent  reference  materials  are  included 
in  four  appendixes.  . 


Preface 


Thia  report  presents  procedures  for  the  systematic  planning  of 
large  river  basin  floodplain  studies.  Work  was  conducted  under  Contract 
No.  CW39-77-M-2001  for  the  U.  S.  Army  Engineer  Waterways  Experiment 
Station  (WES),  Vicksburg,  Miss.,  during  the  period  20  April  1977  through 
30  September  1977.  The  effort  constitutes  part  of  the  Wolf  River  Simu¬ 
lation  Study  being  conducted  for  the  U.  S.  Antty  Engineer  District,  Mem¬ 
phis,  Memphis,  Tenn.,  by  the  WES  under  Appropriation  No.  96X1*902. 

The  3tudy  was  made  and  the  report  prepared  by  Dr.  Gert  Aron,  Depart 
ment  of  Civil  Engineering,  Pennsylvania  State  University,  University 
Park,  Penn.  The  contract  was  monitored  technically  by  Dr.  L.  E.  Link, 
Chief,  Environmental  Research  Sranch  (121E),  Environmental  Systems  Divi¬ 
sion  (ESD),  Mobility  and  Environmental  Systems  Laboratory  (ME3L),  of  WES 
and  Mr.  J.  G.  Collins,  ERB,  Project  Manager,  under  the  general  super¬ 
vision  of  Messrs.  B.  0.  Benn,  Chief,  ESD,  and  W.  G.  Shockley,  Chief, 
MESL. 

The  ESD  is  now  part  of  the  Environmental  Laboratory  of  which 
Dr.  John  Harrison  is  Chief. 

The  Director  of  WES  during  the  period  of  this  study  and  the  prepara 
tion  of  this  report  was  COL  J.  L.  Cannon.  Technical  Director  was 
Mr .  F .  R .  Brown . 
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Conversion  Factors,  U.  S.  Customary  to  Metric  (SI] 
Units  of  Measurement 


U.  S.  customary  units  of  measurement  used  in  this  report  can  be  con¬ 


verted  to  metric  (SI)  units  as  follows: 


Multiply 


cubic  feet  per  second 


inches 

miles  (U.  S.  statute) 
square  miles 


U0U6. 856 
0.02831685 
0.3048 
25.4 

1.609344 

2.589988 


_ To  Obtain _ 

square  metres 

cubic  metres  per  second 

metres 

millimetres 

kilometres 

square  kilometres 
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PLANNING  MANUAL  FOR  SYSTPiATIC 
RIVER  HASIN  FLOODPLAIN  STUDIES 


Introduction 

1.  Floodplain  management  has  become  a  topic  of  high  national 
priority  in  the  last  decade,  and  billions  of  dollars  will  probably  be 
spent  in  the  near  future  on  floodplain  delineations  and  flood  level 
computations. 

2.  Hydrologic  and  hydraulic  process  simulation  computer  programs 
such  as  KEC-1  and  HEC-2  have  greatly  increased  the  efficiency  with 
which  calculations  can  be  performed.*  At  present,  though,  most  flood- 
plain  studies  are  being  conducted  in  parts,  by  first  performing  the  hy¬ 
drologic  analysis,  then  following  up  with  the  hydraulic  computations. 
Much  time  is  spent,  and  possibly  wasted,  in  collecting  data  from  scat¬ 
tered  sources,  running  the  hydrologic  analysis,  then  processing,  the  out¬ 
put  of  this  analysis  into  the  appropriate  input  for  the  hydraulic  com¬ 
putations.  Surveying  and  plotting  of  stream  cross  sections  and 
converting  these  to  the  proper  input  format  for  backwater  programs  is 
another  highly  tedious  Job  subject  to  many  possible  human  input  errors 
or  inconsistencies. 

3.  The  first  objective  of  this  manual  is  to  present  guidelines 
for  standardizing  watershed  subdivision  and  collecting  and  processing 
data  into  unambiguous  formats  that  can  be  converted  automatically  into 
inputs  for  computer  programs  HEC-1  and  HEC-2.  The  second  objective  is 
to  maximize  the  ease  with  which  these  two  programs  con  be  interfaced 
so  that  a  minimum  of  input  data  needs  to  be  introduced  manually  once 
the  main  data  tapes  with  information  on  basin  topography,  soils,  land 
use,  and  stream  cross  sections  have  been  compiled.  A  uniform  treatment 
of  bridge  geometry  will  also  be  suggested,  which,  while  not  necessarily 
using  the  most  appropriate  hydraulic  computational  procedure  in  every 


HEC  *  Hydrologic  Engineering  Center 


case,  should  provide  water  surface  profiles  commensurate  in  precision 
with  presently  used  backwater  computation  procedures. 

Organizational  Step  Sequence 

1*.  A  systematic  hydrologic-hydraulic  floodplain  analysis  should 
involve  the  following  steps: 

a.  Selection  of  Watersheds  to  be  Studied 

b.  Identification  of  Primary  and  Secondary  Areas  in  the 
Watershed 

c_.  Planning  Cross-Section  Location  and  Length 

d.  Subdivision  of  Watershed  and  Subarea  Numbering 

£.  Collection,  Coding,  and  Digitization  of  Cross-Section  and 
Bridge  Data 

f .  Collection  of  Soils,  Topography,  Land  Use,  and  other  Runoff- 
Determining  Data 

£.  Establishment  of  a  Digital  Data  Bank 

h.  Collection  of  Measured  Rainfall-Runoff  Data 

i. .  Model  Calibration 

J_.  Choice  of  Event  Storms 
k.  Generation  of  Event  Hydrographs 
1_.  Backwater  Computations 
a.  Areal  Delineation  of  Floods 
n.  Flood  Damage  Computations 

Selection  of  water- 
sheds  to  be  studied 

5.  The  selection  of  the  watersheds  for  which  floodplain  studies 
should  be  made  depends  primarily  on  the  urgency  of  the  situation.  Obvi¬ 
ously,  a  densely  developed  watershed  with  a  history  of  costly  floods 
merits  preference  over  a  wilderness  area  in  which  flooding  is  an  en¬ 
tirely  acceptable  event. 

6.  There  Is  no  lower  limit  for  a  watershed  size,  although  it  may 
be  too  cumbersome  to  bother  about  an  elaborate  systematic  scheme  for  a 
watershed  of  one  square  mile  or  less.  The  upper  limit  of  a  watershed 
size  to  be  handled  in  a  single  effort  depends  mostly  on  computer 
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storage  capacity  and  available  time  to  complete  the  project. 

7.  Whatever  the  size  of  the  area  to  be  modeled,  it  should  com¬ 
prise  an  entire  watershed.  No  streams  should  cross  the  boundary  except 
at  the  outlet  section.  For  example,  the  method  outlined  should  not  be 
used  to  study  river  levels  between  points  A  and  B  in  Figure  1  unless 
reliable  inflow  hydrographs  for  point  A  have  been  generated  in  a  pre¬ 
vious  study. 


Identification  of  primary  and 
secondary  areas  in  the  watershed 


Given  sufficient  funds,  the  entire  watershed  can  be  considered 


Figure  1.  Partial  watershed  with  separately 
determined  inflow 
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of  primary  importance  and  designated  for  detailed  data  acquisition  and 
modeling.  Even  areas  which  are  complete  widerness  or  swamps  at  present 
may  at  some  future  date  become  attractive  for  development. 

9.  In  practice,  it  will  usually  be  found  that  a  detailed 
hydrologic-hydraulic  modeling  effort  of  an  entire  watershed  of  1000 
square  miles*  or  so  would  become  prohibitively  expensive.  Thus,  some 
areas  may  be  considered  of  secondary  importance  and  given  a  more  super¬ 
ficial  modeling  treatment.  Subareas  may  be  made  larger  than  in  the 
detailed  study  area  and  cross  sections  may  be  spaced  farther  at>art.  It 
i 3  important,  however,  to  leave  provisions  for  later  conversion  of 
these  secondary  areas  to  primary  interest  areas  without  causing  a  major 
disruption  in  the  model  organization.  This  involves  a  flexible  coding 
of  the  subVasins  such  that  future  addition  of  subbasins  will  not  re¬ 
quire  the  sequential  recoding  of  the  entire  system. 

Planning  cross- 

section  location  and  length 

10.  Once  it  is  determined  which  portions  of  the  watershed  are  to 
be  classified  as  primary  and  secondary  interest  areas,  the  locations 
and  lengths  of  the  cross  sections  should  be  planned.  Slight  changes 

in  cross-section  location,  or  the  addition  of  a  few  cross  sections,  can 
always  be  made  at  the  time  of  the  survey. 

11.  First,  all  streams  and  tributaries  of  hydrologic  importance 
must  be  identified.  These  are  the  streams  for  which  event  floods  of 
interest  and  associated  flood  levels  are  to  be  computed.  This  classi¬ 
fication  may  become  a  rather  subjective  task,  strongly  dependent  on 
the  funds  available  for  the  flood  study.  Referring  to  Figure  2,  assume 
that  streams  a,  b,  c,  and  e  have  been  designated  as  hydrologically 
important. 

12.  At  all  Junctions  of  important  streams  two  sets  of  cross- 
section  data  must  be  obtained:  one  section  should  he  located  about  one 
stream  width  below  the  Junction,  another  section  between  one  and  two 


A  table  of  factors  for  converting  U.  S.  customary  units  of  measure¬ 
ment  to  metric  (SI)  can  be  found  on  page  3. 


Figure  2.  Cro3s-section  location  on  a 
small  drainage  basin 


stream  widths  above  the  Junction.  The  upstream  cross  section  must  be 
split  at  the  divide,  as  shown  in  Figure  2  at  the  Junctions  of  streams 
a  and  b,  b  and  c,  and  a  and  e.  This  separation  of  data  into  two  inde¬ 
pendent  upstream  sections  is  absolutely  necessary  for  the  correct  execu¬ 
tion  of  HBC-2  backwater  computations.  The  cross  section  should  extend 
through  the  floodplain  on  both  sides  of  the  stream  along  a  course  nor¬ 
mal  to  the  flow  lines  (i.e.,  approximately  perpendicular  to  the  topo¬ 
graphic  contour  lines),  except  at  a  bridge  crossing  where  the  cross  sec¬ 
tion  may  follow  the  road  bed. 

13.  Cross  sections  are  required  at  all  bridges.  The  number  and 
spacing  of  these  cross  sections  will  be  discussed  in  more  detail  in  the 
bridge  data  instructions,  step  e.  Cross  sections  should  also  be  pro¬ 
vided  at  all  control  locations,  like  weirs  or  other  abrupt  flow  con¬ 
strictions,  at  which  critical  flow  is  expected  to  occur  for  some  or 
all  of  the  discharges  of  interest.  The  same  holds  for  locations  of 
stream  gages  which  might  be  used  for  the  calibration  of  backwater  com¬ 
putations,  as  well  as  those  points  at  which  the  channel  and/or  overbank 


8 


geometry  changes  appreciably.  Control  sections  are  often  not  apparent 
on  the  topographic  maps;  therefore,  it  would  be  extremely  useful  if  a 
hydrologist  or  iiydraulic  engineer  accompanies  the  survey  crew,  possibly 
as  a  crew  member,  or  at  least  performs  a  field  reconuissance  to  locate 
all  cross  sections. 

1U.  Observance  of  the  above  set  of  guidelines  for  cross-section 

1  2 

location  or  a  similar  set  provided  by  HEC  '  will  usually  result  in 
rather  close  spacing  of  cross  sections  in  developed  areas  and  large 
spacings  in  less  developed  and  hydrologiealiy  less  important  areas. 

This  brings  up  the  problem  of  maximum  allowable  reach,  length,  an  econom¬ 
ically  important  subject,  especially  in  relatively  inaccessible  stream 
reaches,  where  cross-section  surveys  are  expensive. 

15.  Barr  Engineering  Co.'  cite3  a  rule  originating  in  Minnesota 
that  makes  the  maximum  reach  length  L  a  function  of  stream  slope  S 
as  follows: 

:'lope,  ft /mile  Maximum  Reach  Length,  ft 

<2  2,61*0 

2-3  1,800 

>3  •  1,200 

This  rule,  when  appl'.ed  to  some  stream  data,  war.  found  overly  conserva¬ 
tive,  especially  In  wide  and  mature  rivers  like  the  Mississippi  where 
reach  lengths  of  2  miles  or  longer  might  be  quite  adequate. 

16.  Tavener  performed  an  analytical  study  of  backwater  computa¬ 
tional  errors  in  natural  stream  channels,  which  resulted  in  an  equation 
for  majcimum  reach  length  as  a  function  of  the  stream,  gradient,  the 
Froude  Somber,  wetted  cross-section  area,  ar.d  wetted  perimeter,  as  well 
as  the  top  width  of  the  stream  at  any  section.  Aron  modified  the 
equation  and  converted  it  to  a  nomograph  solution  which  is  less  tedious 
to  derive.  Both  the  Tavener  equation  and  the  Aron  nomograph,  however, 
require  dimensions  which  are  not  known  until  after  flood  profiles  have 
been  letenalned;  therefore,  both  of  these  methods  are  suited  chiefly  to 
check  whether  a  chosen  cross-section  spacing  was  adequate.  Appendix  D 
presents  formulas  for  backwater  computations. 
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17.  As  part  of  the  research  leading  to  this  manual,  backwater  com¬ 
putations  were  run  on  a  hypothetical  channel  of  parabolic  cross  section 
and  randomly  varying  width  and  slope.  Channel  width  was  varied  between 
20  and  100  ft,  depth  between  1.5  and  U  ft,  slope  between  0.0006  and 
0.02,  and  discharge  between  100  and  1000  cfs.  With  any  combination  of 
these  variables,  reach  length  was  varied  between  POO  and  «800  ft.  Using 
200-ft  reach  length  results  as  a  standard  for  comparison,  water  profile 
errors  occurred  sporadically,  and  no  consistent  formula  for  maximum 
reach  length  could  be  found.  This  was  because  errors  due  to  expansion 
and  contraction  losses  interfered  randomly  with  errors  due  to  friction 
losses.  Examining  only  the  largest  errors  in  computed  water  surface 
elevation,  these  seemed  to  be  related  to  the  irregularities  of  the  chan¬ 
nel  slope.  The  equation  developed  from  the  observations  Is 


max 


"  .  ^ 

S  ♦  A  S 


(1) 


where 

L  *  the  maximum  advisable  reach  length,  ft  or  a 
MI 

e  ■  the  largest  acceptable  error  in  computed  water  surface 
elevation,  ft  or  a 

S  *  the  average  channel  slope,  ft/ft 

AS  *  the  range  of  channel  slopes,  ft/ ft 

in  the  vicinity  of  the  Intended  reach,  as  determined  fro®  a  detailed 
topographic  map  or  a  preliminary  survey  of  low  water  surface  profiles. 
The  streombed  itBelf  may  have  local  scour  holes  and  irregularities 
which  have  negligible  effect  on  the  computed  flood  profiles. 

18.  For  a  sample  application,  take  a  stream  with  an  average  chan¬ 
nel  slope  of  2  ft/mile  over  the  lover  5  miles.  According  to  the  Barr 
rule,  cross  sections  should  be  taken  at  about  one  half  mile  intervals. 
If  the  maximum  allowable  error  was  0.1*  ft  and  the  slope  varied  between 
1  and  3  ft/mile.  Equation  1  would  give  the  same  maximum  reach  length  of 
one  half  mile.  If,  on  the  other  hand,  a  1.0- ft  error  was  acceptable 
and  the  slope  varied  only  between  1.5  and  2.5  ft/oile,  reach  lengths  as 
large  as  1.7  miles  would  be  acceptable.  The  use  of  Equation  1  must. 
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however,  be  subject  to  the  previously  stated  rules  for  cross-section 
location.  It  should  also  be  kept  in  mind  that  the  computations  leading 
to  Equation  1  vere  made  exclusively  for  subcritical  flow  conditions. 

19.  Under  certain  conditions,  reach  lengths  used  may,  at  least 
partially,  depend  on  the  friction  slope  average  method  selected;  thiB 
topic  is  discussed  herein.  A  number  of  equations  have  been  proposed  for 
calculating  friction  losses,  four  of  which  are  options  in  the  HEC-2 
computer  program.  These  four  equations  are  listed  below: 

Average  Conveyance  Equation 


(2a) 


Average  Eric*:  .'Icpe  Equation 


HL  -  L 


S,  ♦  S„ 
1  £ 


•cocetri  c  Mean  Friction  .'lope  Equation 

KL  •  L  (B^g)0,5 


(2b) 


(2c) 


Harmonic  V'-w  in  .  -pc  Equation 


HL  -  L 


2  3,  S_ 


(2d) 


where 

HL 


s 

3, 


friction  loss  for  reach 
discharge-weighted  reach  length 
total  discharge  at  upstream  end  of  reach 
total  discharge  at  downstream  end  of  reach 
total  conveyance  at  upstream,  end  of  reach 
total  conveyance  at  downstream  end  of  reach 
friction  slope  at  upstream,  end  of  reach 
friction  slope  at  downstream  end  of  reach 
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20.  Equation  2a  is  used  in  all  KEC-2  source  decks  dated  August 
1971  that  contain  Modification  56.  It  is  presently  recommended  by  KEC 
for  all  flov  conditions  until  research  shovs  that  another  equation  is 
clearly  more  applicable. 

21.  Equation  2b  represents  the  arithmetic  average  of  the  friction 

4 

slopes  S  and  S-  and  has  been  recommended  by  Tavener  as  the  prefer 

12  g 

able  equation  for  M-l  curves.  Reed  and  Volfkill  confirmed  Tavener's 
recommendations . 

22.  Equation  2c  represents  the  geometric  average  of  the  friction 

slopes  and  S„  and  is  presently  used  in  the  U.  S.  Geological 

Survey  (USGS)  computer  program  for  calculating  profiles. 

23.  Equation  2d  is  equivalent  to  averaging  the  inverse  values  of 

slopes  and  S„  and  has  been  recommended  by  Tavener  for  M-2 

profiles. 

24.  Modification  50  of  the  computer  program  HEC-2,  dated  November 
1976,  allows  the  choice  among  Equations  2a  to  2d  by  specifying  a  vari¬ 
able  1HLBQ  in  field  1  of  input  card  J 6,  as  listed  below: 


Value  of 
IHLEQ 

0 

1 


Friction  Loss  Equation 


Equation  ?a  is  used 

Program  selects  equation  based  on  flow 
conditions 


2  Equation  2b  in  used 

3  Equation  2c  is  used 

4  Equation  2d  is  used 

If  IHLEQ  is  set  equal  to  1,  the  program  selects  a  friction  loss  equa¬ 
tion  for  a  reach  in  accordance  with  criteria  below: 


PrafUs,  TOrre 


Subcrltical  (Ml, SI) 
Subcritical  (M2) 
Supercritical  (S2) 
Supercritical  (M3, S3) 


Is  Friction  Slope  at  Current 
Cross  Section  Greater  than 
Friction  Slope  at  Preceding  Equation 
_ Crsaa  SccUoaS _  — 


Yes 

2b 

No 

2d 

YC8 

2b 

No 

2c 

12 


2 5.  At  present,  HEC’s  stand  on  friction  slope  averaging  is  that 
the  choice  of  a  particular  equation  nukes  little  difference  in  the  com¬ 
puted  water  surface  profiles,  as  quoted  below. 

Experience  to  date  indicates  that,  application  of 
criteria  in  Table  1  produces  water  surface  profiles 
that  only  rarely  differ  by  more  than  0.2  feet  fro® 
profiles  determined  with  Equation  (2a).  In  a  few  in¬ 
stances,  application  of  criteria  in  Table  1  enabled 
determination  of  'balanced'  water  surface  elevations 
at  cross  sections  for  which  Equation  (2a)  could  not 
produce  a  solution  to  the  energy  equation.  Any  of 
the  alternative  friction  loss  equations  will  produce 
satisfactory  estimates  provided  that  reach  lengths 
are  not  too  long.  The  advantage  that  is  sought  in 
alternative  friction  loss  formulations  is  to  be  able 
to  maximize  reach  lengths  without  sacrificing  profile 
accuracy. 

Two  specific  situations,  however,  are  Illustrated  in  Figures  3  and  k  in 
which  the  improper  choice  of  a  slop>e-ave raging  equation  can  lead  to 
major  errors  in  computed  water  surface  elevation. 

26.  For  the  arithmetic  slope-averaging  method  (Fquation  2b),  the 
most  critical  condition  is  illustrated  in  Figure  3;  here  the  downstream 
slope  is  abnormally  high,  due  perhaps  to  a  combination  of  a  narrow 
floodway  and  high  channel  roughness.  If  in  this  typical  M2  profile 

a  few  short  reaches  are  used  Just  upstream  of  the  construction,  the 
steep  energy  slope  is  gradually  reduced  to  a  value  which  approximates 
normal  flow  conditions  in  the  channel.  With  a  long  reach,  on  the  other 
hand,  ‘he  backwater  iterations  may  converge  to  a  very  large  upstream 
depth  because  the  minimum  average  friction  slope  used  for  the  reach 
cannot  be  smaller  than  one  half  of  the  downstream  friction  slope,  which 
can  be  as  high  as  0.02  at  or  near  critical  flow  conditions  in  shallow 
streams . 


27.  Equation  2a,  presently  preferred  by  HEC,  produces  adequately 
accurate  flood  profiles  in  the  M2  profile  situation  described  above 
even  with  reaches  of  several  hundred  feet  in  length,  but  can  cause 
problems  in  the  Ml  profile  situation  shown  in  Figure  fc,  in  which  the 
flood  water,  backed  up  by  a  bridge  or  other  obstruction,  occupies  a 
wide  and  deep  downstream  section  resulting  in  a  very  large  conveyance. 
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Computed  Water  Surface 


Realistic  Water  Surface 
with  M2  Profile  — ' 


Water  profile  error  due  to  use  of  incorrect  slope 
averaging  equation  used  on  M2  profile 

Realistic  Water  Surface  with 
Ml  Profile 


Figure  3 


Computed  Water  Surface 


Bridge  Deck 


Water  profile  error  due  to  use  of  incorrect  slope 
averaging  equation  used  on  Ml  profile 


This  leads  to  a  condition  under  which,  no  matter  how  narrow  the  upstream 
section  is,  the  average  conveyance  cannot  be  lower  than  one  half  the 
downstream  conveyance.  The  consequence  of  such  a  situation  is  an  in¬ 
escapably  low  friction  slope  which  forces  the  HEC-2  program  Into  a 
critical  depth  default  if  the  reach  length  is  too  long. 

28.  Because  of  the  observations  made  in  paragraphs  26  and  27,  it 
is  recommended  that  input  card  J6  be  used  with  variable  IHLBQ  set  equal 
to  1. 


Subdivision  of  water¬ 
shed  and  subarea  number inj 


29.  The  segmentation  or  subdivision  of  a  watershed  is  a  eoranon 
feature  of  all  major  storm  water  runoff  models.  Ideally,  it  would  be 


convenient  to  have  all  watershed  subareas  of  equal  size,  but,  as  shown 
in  Figure  5,  the  location  and  density  of  subdivision  boundaries  are 
primarily  dependent  on  the  density  of  streams  important  enough  for 
flood  level  computations. 

30.  In  the  watershed  subdivision  process,  a  hierarchal  system  may 
be  used.  The  upper  hierarchy  is  obviously  the  watershed  to  be  modeled. 


Figure  5.  Two-level  watershed  subdivision  into 
sufcbasins  and  minibasins,  using  numbering  se¬ 
quence  alternative  No.  1 


The  watershed  may  be  divided  into  subbasins  and  these  further  into  mini¬ 
basins.  In  small  to  medium-size  watersheds  in  which  the  number  of  sub¬ 
divisions  is  sufficiently  small  to  be  processed  in  one  single  HEC-1  run, 
the  watershed  may  be  divided  directly  into  minibasins  in  what  will  be 
defined  hereafter  as  a  one-level  subdivision  process.  The  convenience 
of  dividing  a  larger  watershed  in  a  two-level  process  into  Bubbasins 
and  minibasins  will  be  discussed  below  in  connection  with  the  subarea 
number  coding.  The  term  subarea  will  be  used  in  describing  properties 
or  processes  common  to  subbasins  and  minib&sins. 

31.  The  first  rule  in  subdividing  a  watershed  must  be  that  no  two 
streams  of  hydrologic  importance  can  be  located  in  the  same  minibasin. 

A  logical  procedure  in  watershed  subdivision  is  to  follow  the  main 
stream  of  the  watershed,  starting  at  the  outlet  section  and  proceeding 
upstream,  outlining  the  watersheds  of  all  the  tributaries  of  hydrologic 
importance  as  defined  previously.  This  involves  drawing  boundary  lines 
along  the  divides  between  adjacent  streams  of  hydrologic  importance. 

If  a  subarea  contains  a  reservoir  large  enough  to  be  included  in  the 
HEC-1  modeling  operation,  its  drainage  area  should  also  be  separated 
as  a  subarea.  Typical  minibasin  boundary  lines  are  3hcwr.  in  Figure  6. 

32.  As  a  general  rule,  a  minibasin  should  not  be  larger  than  one 
square  mile  in  densely  developed  &rea3,  but  in  vilderness  areas  which 
are  highly  uniform  or  for  which  little  hydrologic  information  is  avail¬ 
able,  subareas  can  be  as  large  as  10  square  miles.  If  the  Soil  Con¬ 
servation  Service  (SCS)  unit  hydrograph  generating  procedure  is  to  be 
used,  however,  minibasins  should  not  be  excessively  large,  perhaps  less 
than  five  square  miles,  because  the  SCS  has  tested  their  procedure  only 
on  areas  smaller  than  2000  acres.  It  is  also  advisable  to  subdivide 
areas  of  highly  nonsymetric  shape  for  which  travel  times  of  different 
portions  would  otherwise  vary  greatly.  When  a  subbasin  must  be  sub¬ 
divided,  divisions  should  be  made  at  surveyed  cross-section  locations 
if  possible.  A  bridge  is  a  good  location  for  a  subdivision  line,  for 
example. 

33.  In  labeling  the  subareas  and  cross  sections,  a  numerical 
rather  than  alphabetic  system  should  be  used,  because  computer  programs 
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Figure  6.  One-level  watershed  subdivision 
into  minibasins 


can  sort  numerical  information  with  much  more  ease.  Subareas  may  be 
numbered  in  an  upstream  or  downstream  sequence.  The  upstream  numbering 
sequence  has  the  advantage  that  the  downstream  limit  of  a  drainage  basin 
is  usually  the  point  of  outflow  into  the  receiving  water  body  and  thus 
is  well  defined.  At  the  upstream  end,  on  the  other  hand,  a  stream  usu¬ 
ally  fades  into  a  creek  and  finally  into  a  gully,  and  the  upstream  ex¬ 
tremity  of  a  stream  is  thus  often  difficult  to  define.  During  the 
modeling  effort  the  hydrologist  may  decide  to  add  or  delete  one  or  more 
upstream  sections  or  basins;  if  so,  fewer  corrections  in  the  numbering 
would  be  necessary  if  an  upstream  sequence  is  used. 

31*.  The  numbering  sequence  must  follow  a  logic  which  can  be  under¬ 
stood  by  a  computer  program.  In  one-level  subdivision  processes  a 
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single  sequence  of  numbers  may  be  used  as  illustrated  in  Figure  6.  The 
largest  of  the  streams  is  chosen  us  the  main  stem,  and  the  minibasin 
numbering  follows  the  sequence  1  to  Wo  the  upper  minibasin  of  this 
stream.  Subsequently,  the  stream  is  scanned  for  tributaries,  starting 
at  the  upper  minibasin,  namely  b,  and  moving  downstream.  At  point  £ 
the  first  tributary  is  found,  consisting  of  only  one  minibasin,  which 
is  labeled  5.  Continuing  farther  down  to  point  b,  a  tributary  flowing 
through  several  minibasins  is  found  and  the  numbering  process  is  re¬ 
peated,  numbering  the  minibasins  6  to  8,  then  scanning  downstream  to 
pick  up  one  more  tributary  at  point  £,  which  requires  only  one  mini¬ 
basin,  numbered  9* 

35.  The  above  described  numbering  sequence  is  compatible  with  the 
computational  sequence  of  HEC-2  and  the  input  requirements  of  HYDPAR2, 

a  program  described  in  Appendix  B,  which  converts  hydrologic,  topo¬ 
graphic,  and  land  use  data  into  an  HEC-1  input  deck. 

36.  In  modeling  larger  watersheds,  it  may  be  advisable  to  divide 
the  watershed  into  subbasins  and  minibasins.  The  reasons  for  using  a 
two-level  subdivision  may  be  the  limit  on  the  number  of  subareas  a 
program  can  handle  in  one  run,  or  the  loss  in  subdivision  flexibility 
due  to  a  large  number  of  sequentially  numbered  minibasins.  With  refer¬ 
ence  to  the  latter,  assume,  for  example,  that  the  watershed  shown  in 
Figure  6  was  divided  into  hundreds  of  minibasins  and  that  the  uppermost 
minibasin  of  the  main  stem  was  labeled  So.  19.  After  numbering  and 
recording  all  other  subbasins  sequentially  it  was  found  necessary  to 
divide  minibasin  19  into  two  parts  which  should  logically  be  renumbered 
19  and  20,  respectively.  This  small  change  then  starts  a  cascade  of 
minibasin  number  adjustments  which  can  become  a  rather  tedious  task  and 
lead  to  con'Usion  between  old  and  new  numbering  sequences.  Figures  5 
and  T  illustrate  cases  of  two-level  subdivision  and  numbering  by  sub¬ 
basins  and  minibasins. 

;I7.  The  subbasin  numbering  may  be  made  to  follow  the  same  logic 
as  prescribed  in  Figure  6  for  the  one-level  subdivision.  This  results 
in  a  numbering  sequence  as  shown  in  Figure  5,  following  the  main  stream 
stem  upward  to  the  uppermost  subbasin,  then  numbering  all  drainage  areas 
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of  the  tributaries  in  downstream  sequence.  17118  numbering  sequence  for 
subbasins  which  are  further  subdivided  looks  awkward;  therefore,  an  al¬ 
ternative  sequence  is  suggested.  In  Figure  7,  drainage  areas  of  tribu¬ 
tary  streams  and  areas  draining  directly  Into  the  main  stem  receive 
equal  subarea  status  and  are  numbered  consecutively  in  upstream  order. 

No  tributary  will  occupy  more  than  one  subbasin,  but  any  subbasin  can  be 
subdivided  into  minibasins. 

38.  The  minibasin  numbering  follows  the  rules  laid  out  in  Figure  6 


Figure  7.  ?wo-level  watershed  subdivision  into 
subbasins  and  minibasins,  using  numbering  se¬ 
quence  alternative  No.  2 


for  single-level  subdivision.  In  this  manner  each  subbasin  can  be 
processed  individually  or  as  part  of  a  group  by  program  HYDPAR2  into  an 
input  deck  for  the  hydrologic  program  HEC-1.  Four  digits  are  reserved 
for  minibasin  numbers;  the  first  two  of  which  identify  the  subbasin,  if 
any,  the  last  two  digits  identify  the  minibasin  within  it. 

39-  Finally  the  cross  sections  within  the  minibasins  must  be  num¬ 
bered.  This  number  should  be  the  section  number  SECNO  entered  in-  field 
1  of  the  XI  card:;  of  the  HEC-2  input  deck.  SECNO  is  restricted  to  3ix 
digits,  the  first  of  which  may  be  needed  for  a  minus  sign.  Therefore, 
the  section  numbers  consist  of  the  four  digits  identifying  the  mini¬ 
basin,  followed  by  one  digit  which  identifies  the  cross  sections  within 
the  minibasin.  There  will  rarely  be  more  than  10  cross  sections  in  a 
mlnlbaain,  except  where  several  bridges  cross  the  stream  in  close  prox¬ 
imity,  perhaps  requiring  a  total  number  of  regular  plus  auxiliary  sec¬ 
tions  larger  than  10. 

U0.  Figure  8  is  an  enlargement  of  the  lower  portion  of  subba3in  2 
in  Figure  7.  It  war  found  convenient  to  label  any  cross  section  imme¬ 
diately  above  a  downstream  minibasin  boundary  with  the  digit  "0"  and 
those  immediately  below  the  upstream  boundary  with  the  digit  "9";  be¬ 
tween  these  boundaries  any  digits  between  1  and  8  may  be  assigned, 
which  provides  the  flexibility  of  later  insertion  of  additional  cross 
sections.  Neither  the  program  HEC-2  nor  XSEC-PHOC,*  which  generates 
HEC-2  input  decks  from  cross-sectional  data,  requires  an  unbroken  se¬ 
quence  of  cross  section  numbers;  HYDPAP2,  however,  recognizer,  only  a 
section  number  ending  with  9  as  legitimate  point  at  which  the  flows 
from  two  upstream  tributaries  can  be  combined. 

Ul.  Auxiliary  sections  occur  at  bridges  and  other  locations  at 
which  a  surveyed  section  is  repeated  upstream  or  downstream.  These 
auxiliary  sections  are  labeled  using  the  surveyed  section  number  with¬ 
out  the  subbasin  number,  followed  by  a  decimal  point  and  decimals  1, 
2,...  in  upstream  sequence.  As  an  example,  the  second  auxiliary  sec¬ 
tion  at  a  Section  Ho.  3218U  could  be  numbered  18U.2. 


•  See  Appendix  A. 
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Figure  8.  Cross  section  numbering  within  minibasins 


Collection,  coding,  and 
digitization  of  cross- 
aection  and  bridge  data 


r 
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a.  Location  of  cross  section,  generally  measured  as  the  dis¬ 
tance  along  the  stream  channel  from  the  stream's  outflow 
point.  This  distance  may  be  measured  in  feet,  miles,  or 
kilometres. 

b.  Angle  by  which  a  cross  section  deviates  from  the  normal 
to  the  stream  channel,  expressed  in  degrees.  If  this 
angle  is  less  than  20  deg,  it  may  be  ignored. 

£.  Distance  by  which  the  left  and  right  overbank  reach 
lengths  between  sections  differ  from  the  reach  length 
measured  along  the  channel.  May  be  measured  from  maps 
and  expressed  in  feet  or  metres. 

Number  and  location  of  auxiliary  cross  sections  created 
by  transferring  surveyed  cross  sections  a  stated  distance 
upstream  (♦)  or  downstream  (-). 

£.  Discharges  as  they  vary  with  distance  from  the  stream  out¬ 
let.  These  estimates  usually  come  from  separate  hydro- 
logic  analyses  and  are  expressed  in  cubic  feet  or  cubic 
metres  per  second.  A  set  of  1  to  lL  design  discharges  is 
specifies  ahead  of  the  first  cross  section,  and  these  dis¬ 
charges  will  remain  constant  until  superceded  by  a  new 
set  of  discharges  entered. 

£.  Location  and  height  of  real  or  artificial  levees.  Whether 
levee  instructions  should  be  generated  for  HEC-2  depends 
on  the  ease  with  which  the  water  can  flow  into  and  out  of 
low-lying  overbank  areas.  (See  HK-2  Users  Manual  for 
levee  card  instructions. ) 

£.  Left  and  right  main  channel  bank  stations.  These  are  the 
limits  within  which  the  streamflow  is  computed  with  the 
ln-bank  roughness  factor. 

h.  Bridge  width  and  height  of  railings.  These  data  are  fre¬ 
quently  omitted  in  cross-section  data.  The  height  of  the 
railings  should  be  expressed  as  the  solid  height.  For 
example,  if  the  railings  consist  of  three  L-in.  horizon¬ 
tal  pipes,  the  height  of  the  railing  should  be  given  as 
1.0  ft. 

i_.  Stations  along  a  bridge  section  within  which  weir  flow 

is  expected  to  occur  when  the  water  surface  elevation  ex¬ 
ceeds  that  of  the  bridge  deck.  These  are  needed  only 
when  the  Special  Bridge  routine  is  to  be  U3ed.  Sometimes, 
a  bridge  cross  section  may  be  a  mile  long.  Obstructions 
might  well  prevent  the  flow  from  freely  moving  to  a  point 
one  half  mile  from  the  stream  and  cascading  across  the 
road.  The  engineer  must  make  a  Judgment  decision  regard¬ 
ing  the  total  width  of  the  weir  flow.  This  width  may  be 
expressed  as  the  first  and  last  stations  for  which  BT 
cards  should  be  created,  or  simply  as  a  width,  in  feet  or 
metres,  to  be  entered  in  field  b  of  the  SB  card  of  HEC-2. 
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J_.  Manning's  roughness  factors  n,  as  well  as  flow  contraction 
and  expansion  coefficients.  Wherever  a  change  in  n-values 
i s  deemed  desirable,  either  three  such  values,  for  left 
overbanX,  right  overbanX,  and  channel  may  be  listed,  or 
more  than  three  values  and  their  respective  end  stations. 

Uh.  Options  which  can  be  used  include  the  choice  of  having  all 
input  data  printed  out,  of  adding  extra  comment  lines  to  the  input,  and 
of  using  feet,  miles,  or  Xilometrea  as  a  measure  of  the  sections'  dis¬ 
tance  fro®  the  watershed  outlet. 

1.  .  A  computer  program.,  XS1T-PR0C,  wn  written  to  convert  the  con¬ 
tents  of  a  cros3-section  tap>e  file  into  standard  HEC-2  input  cards.  The 
processes  followed  by  the  program  are  shown  in  Figure  A1 .  Detailed  in¬ 
put  instructions  to  XSEC-PROC,  as  well  as  a  program,  statement  listing, 
a  list  of  definitions  of  program  variables,  and  a  sample  problem  with 
input  and  output  lines,  are  included  in  Appendix  A. 

Collection  of  soils, 
topography,  land  use,  and 
other  runoff-determining  data 

!»6.  Physical  watershed  and  subarea  data  liXe  land  slopes,  stream 
lengths,  surface  roughness  factor:,  and  soil  properties  affecting  run¬ 
off  are  needed  for  the  development  of  design  hydrographs. 

•*7.  Quadrangle  maps  by  USC.S,  SCS  soil:;  and  lar.i  use  maps ,  and 
aerial  and  photographs  or  inftared  imagery  are  prime  sourcer  of  basic 
inf  rmatlon.  for  modeling  snal 1  watersheds,  a  manual  procedure  of  ex¬ 
tracting  watershed  properties  data  and  converting  them  into  HEC-1  input 
is  not  a  difficult  tasX,  but  for  larger  watersheds  a  acre  automated  pro¬ 
cess  of  producing  these  rather  voluminous  KEC-1  input  decXs  is  needed. 
The  input  card  sequence  is  highly  repetitive  and  conducive  to  human 
error  if  produced  in  large  quantities,  but  it  follows  a  simple  and  con¬ 
sistent  logic  amenable  to  computer  processing. 

i*8.  Automatic  interpretation  of  remote  sensing  imagery  has  made 
considerable  progress  over  the  past  two  decodes.  The  U.  S.  Army  Bigi- 
neer  Waterways  Experiment  Station  (WES),  S'C>,  as  well  as  many  univer¬ 
sities,  have  developed  techniques  to  scan  imagery  and  classify  the  con¬ 
tents  by  land  and  topographic  parameters. 
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1*9.  A  program  HYDPAE  has  been  developed  by  the  HEC  for  the  book¬ 
keeping  of  physical  data  and  for  producing  representative  average  values 
of  HEC-1  input  parameters.  The  program  was  expanded  by  WES  into  the 
version  HYDPAR2,  for  which  a  description  and  a  set  of  input  instructions 
are  given  in  Appendix  B.  This  program  scans  a  digital  data  bank  for 
subarea  boundaries  and  watershed  parameters  on  a  grid-by-grid  basis, 
sorts  the  information  by  subbasin  or  minibasin,  and  produces  an  almost 
complete  HEC-1  input  deck  for  a  watershed  or  subbasin. 

Establishment  of 
a  digital  data  bank 

>0.  A  digital  data  bank  for  storing  all  quant itative  data  on 
watershed  geometry,  soil,  and  land  ue  is  indi-.pensibic  in  any  auto¬ 
matic  hydrologic  and  hydraulic  modeling  process. 

51.  One  or  more  standard  geometric  grid  sizer  should  be  chosen 
for  the  data  bank.  In  presently  ongoing  research  at  WES,  a  200-  by 
200-a  grid  network  was  established  for  general  watershed  data,  and  a 
50-  by  t>0-«  network  was  established  for  detailed  land  use  data  vithin 
flood-prone  areas. 

52.  A  standard  geographic  coordinate  system  should  be  eh.  sen  and 
available  caps  rectified  to  agree  with  it.  In  its  mapping,  WES  has 
adopted  the  Universal  Transverse  Mercator  (UTM)  system.  Air  photos  and 
other  imagery  used  must  also  be  rectified  to  agree  with  this  coordinate 
system  so  that  data  from  maps  and  Imagery  can  be  stored  in  the  proper 
grid  spaces .  Subarea  boundaries  are  digitized  and  associated  identifi¬ 
cation  numbers  specified. 

53.  A  program  must  be  written  to  scan  through  the  grid  network, 
identify  each  grid  space  by  the  subbacin  or  ninibasin  to  which  it  be¬ 
longs,  and  reproduce,  upon  request,  a  map  of  all  areas  and  r ubareas  with 
corresponding  environmental  information.  These  maps  should  then  be 
verified  by  visual  comparison  with  the  original  maps  and  other  imagery. 

51*.  After  verification,  the  data  bank  is  ready  to  be  accessed  by 
the  program  HYPPAR2  which  combines  the  watershed  data  and  some  sepa¬ 
rately  dictated  routing  instructions  to  generate  all  but  the  introduc¬ 
tory  input  cards  for  HEC-1. 


Collection  of 

measured  rainfall -runoff  data 


55-  The  final  objective  of  a  hydrologic  and  hydraulic  flood  simu¬ 
lation  model  is  to  reproduce  the  watershed  response  to  storms  of  varying 
recurrence  frequency  under  present  and  future  conditions.  To  verify 
the  model  accuracy,  however,  reliable  data  on  recorded  historic  3tcrm 
and  runoff  events  should  be  used. 

56.  The  nearest  rain  gages  with  reliable  records  should  be  lo¬ 
cated.  The  hyetographs  to  be  used  in  the  model  calibration  can  then  be 
weighted  by  some  appropriate  method  to  reflect  the  extent  to  which  each 
rain  gage  represents  rainfall  in  any  subarea.  Three  available  rainfall¬ 
weighting  procedure:  are  Isohyetal  Contour  Mapping,  Thiesscn  Diagraming , 
and  Weighting  by  the  Inverse  of  the  Distances  to  all  adjacent  gages. 
These  methods  are  applied  to  the  watershed  shown  in  Figure  9.  The 
Isohyetal  Mapping  technique  in,  however,  applicable  only  to  total  storm: 
and  not  to  hyetograph3.  Besides,  thin  method  is  difficult  to  perform 
in  any  computerised  form. 


Figure  9.  Watershed  with  three  rain  gages  and 
Isohyetal  and  Thiessen  Lines 


57.  The  Thiessen  Diagram  weights  rainfall  at  any  point  exclu¬ 
sively  in  favor  of  the  nearest  gage,  even  if  the  difference  between  dis¬ 
tances  to  nearby  rain  gages  is  only  a  few  feet.  This  can  lead  to  abrupt 
changes  in  weighted  rainfall  between  two  adjacent  subareas.  If  the 
watershed  is  divided  into  a  Tine  grid  system  for  land  use  mupping  as 
described  in  paragraph  51,  the  hyetograph  recorded  at  the  rain  gage 
nearest  to  each  grid  space  can  be  assigned  to  that  space  and  the  re¬ 
sulting  grid  space  hyetograph3  averaged  over  the  respective  subarea. 

58.  A  rainfall-weighting  procedure  recently  proposed  and  ideally 
adapted  to  computer  processing  is  that  of  weighting  by  inverse  dis¬ 
tances.  Ac  an  example,  let  point  F  in  Figure  9  represent  the  center  of 
a  grid  space  registered  in  the  digital  data  base.  The  weighted  average 
rainfall  R  amount  for  any  given  time  interval  is 


R  - 


(3) 


where 

P.  *  rainfall  recorded  at  gage  1 

Tj  ■  distance  between  gage  i  and  point  P 
n  *  the  number  of  gages  affecting  average  rainfall  at  point  P 
k  *  exponent 

g 

59-  Dean  and  Synder  conclude  from  the  results  of  recent  research 
that  the  number  n  of  nearest  gages  to  any  point  should  be  three  in 
general,  and  four  in  the  case  of  a  very  dense  rainfall  network.  Dean 
also  recommends  a  value  of  k  8  2  for  the  exponent  in  Equation  3, 
which  in  his  tests  gave  a  slightly  greater  accuracy  in  representing 
rainfall  distributions  and  reduced  the  discrepancy  between  weighted 
rainfall  obtained  by  using  the  three  or  four  rain  gages  nearest  to  a 
point.  A  program  RAINVF.IGHT  was  written  at  the  Department  of  Civil 
frigineering  of  Pennsylvania  State  University,  which  compute:  weighted 
average  rainfall  amount  for  a  small  grid  system  by  the  inverse  distance 
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formula.  The  program  description,  listing,  and  input  instructions  are 
presented  in  Appendix  C. 

60.  A  potentially  serious  drawback  of  any  weighting  method  applied 
to  the  ordinates  of  iiyetcgraphs  has  come  to  light  during  modeling  ef¬ 
forts  performed  by  WES  on  the  Wolf  River  watershed  in  Tennessee.  Unless 
the  hyetographs  recorded  at  the  various  rain  gages  are  synchronized  and 
of  similar  shape,  weighted  averaging  cf  several  hyetographs  may  lead  to 
a  marked  attenuation  of  rainfall  intensities.  Figure  10  shows  hypo¬ 
thetical  hyetographs  from  three  stations  weighted  by  30,  20,  and  50 
percent,  respectively.  Each  of  the  three  hyetographs  has  a  distinct 
peak  describing  a  storm  which  seems  to  be  traveling  from  gage  1  toward 
gage  3.  The  weighted  average  hyetograph  on  line  d  combines  all  peaks 
into  a  single,  drawn-out  storm  of  low-to-nedium  intensity.  This  should, 
at  least  locally,  result  in  highly  reduced  flood  peak  estimates. 

61.  To  maintain  the  peakedness  of  the  recorded  hyetographs  and 
still  give  an  appropriate  weight  to  several  nearby  rain  gages,  the  fol¬ 
lowing  procedure  is  suggested.  For  each  minibasin  to  be  modeled,  the 
three  gages  nearest  to  the  basin',  center  of  gravity  are  found,  and 
the  three  hyetographs'  total  rainfall  amounts  as  well  as  their  temporal 
centers  of  gravity  are  computed.  The  basin  is  then  assigned  a  hyeto¬ 
graph  shape  like  that  of  the  nearest  rain  gage.  Weighting  factors  de¬ 
termined  by  the  Inverse  Distance  or  Thiessen  method  are  then  applied 

to  adjust  rainfall  ordinates  and  timing.  Line  e  in  Figure  10  demon¬ 
strates  this  procedure;  the  peakedness  of  the  storm  recorded  at  gage  3 
is  retained,  but  its  ordinates  are  reduced  and  its  timing  advanced 
according  to  the  weights  assigned  to  gages  1,  2,  and  3.  If  in  addition 
to  some  recording  rain  gages,  nonrecording  gages  are  found  in  the  vicin¬ 
ity  of  the  minibasin,  the  total  storm  amount  from  these  gages  may  be 
used  to  further  weight  the  hyetograph  ordinates  without  changing  its 
shape  and  timing. 

62.  For  needed  calibration  of  the  hydrologic  model,  several  re¬ 
corded  events  of  rainfall  and  runoff  should  be  chosen.  If  possible, 
isolated  Intense  sterms  of  short  duration  should  be  used  because  these 
lend  themselves  best  to  unit  hydrograph  development. 


Rain  gage  1 


Rain  gage  2 
Weight  -  0.2 


Rain  gage  3 


<i.  Conventionally 

weighted  hyetograph 


e.  Hyetograph  weighted 
by  pattern  conserving 
method 


Time  Units 


Figure  10.  Example  of  hyetograph  weighting  methods 


■OT 


63.  A  major  problem  in  the  collection  of  hydrologic  data  is  usu¬ 
ally  the  lack  of  streamflow  records  on  a  watershed.  Ideally,  each  sub¬ 
basin  in  a  watershed  should  have  streamflow  records  in  order  to  cali¬ 
brate  the  rain fall -runoff  model  with  any  degree  of  confidence;  in 
reality  it  is  often  found  that  no  stream  gage  record  whatsoever  is 
available  for  an  entire  watershed. 

Model  interfacing  and  calibration 

6U.  Once  the  watershed  subdivision  has  been  accomplished  and  an 
adequate  amount  of  rainfall-runoff  data  collected,  calibration  modeling 
can  begin.  On  watersheds  with  narrow  valleys  in  which  floodplain 
storage  has  a  small  attenuation  effect  on  flood  hydrographs,  the  model¬ 
ing  process  can  proceed  along  a  simple  sequence  of  operational  steps 
without  any  feedback  loops,  as  shown  in  Figure  11.  The  watershed  data 
i s  processed  by  HYBPAR  and  combined  with  rainfall  data  to  compile  the 
input  deck  for  program  REC-1.  The  Muskingum  routine  may  be  used  for 
hydrograph  routing  with  the  proper  coefficients  determined  from  hydro- 
logic  handbooks.  The  peak  flows  of  the  hydrographn  generated  at  the 
outflow  section  of  each  minibasin  are  then  extracted  for  input  at  the 
appropriate  locations  of  the  HEC-2  input  deck,  which  has  been  largely 
produced  by  program  XSFC-PROC.  Finally,  the  HEC-2  program  is  run  to 
compute  flood  levels;  computed  flows  as  well  as  flood  profiles  may  then 
be  compered  with  observed  flows  and  flood  levels. 

65.  In  a  watershed  with  wide  floodplains,  storage  is  an  important 
factor  in  attenuating  hydrographs.  This  storage,  however,  cannot  be 
known  prior  to  flood  profile  computations.  The  HEC-2  program  may  be 
used  to  compute  cumulative  floodplain  storage  and  offers  the  option  of 
punching  out  storage-outflow  cards  for  HEC-1 ,  but  HEC-2  depends  on  the 
HEC-1  runs  and  output  for  flow  values.  In  this  case  the  looped  process 
sequence  shown  in  Figure  12  should  be  used.  In  this  sequence,  a  pre¬ 
liminary  run  of  HYDFAR2  is  performed  to  generate  a  set  of  synthetic 
streamflow  distributions  based  on  the  equation 

Q  -  c  A°'7  (I*) 
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Figure  11.  Flov  chart  of  systematic  floodplain  and  damage  study 
(narrow  floodplains  with  small  storage  effects) 


where 

Q  *  peak  flow  at  a  given  point  on  the  watershed 

A  *  drainage  area  upstream  from  this  point 

c  B  a  constant  given  10  values  between  zero  and  a  value  resulting 
in  the  maximum  flood  possibly  expected  at  the  watershed  outlet 
as  determined  by  some  chosen  conservative  empirical  formula  or 
available  flood  frequency  data. 

66.  Program  HE C-2  is  then  run  with  these  synthetic  flows  to  pro¬ 
duce  a  table  of  storage-outflow  values  for  HEC-1.  Subsequently  HYPPAR2 
is  called  upon  to  generate  the  HEC-1  input  deck,  HEC-1  is  run,  and  its 
output  used  to  make  the  final  run  of  HEC-2. 

67.  In  calibrating  the  models,  the  hydrographs  generated  by  HEC-1 
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Figure  12.  Flow  chart  of  systematic  floodplain  and  damage  study 
(broad  floodplains  with  major  storage  effects) 


should  be  compared  with  recorded  hydrographe  wherever  the  latter  are 
available.  Model  parameters  may  be  adjusted  as  follows:  if  the  com¬ 
puted  hydrograph  volumes  are  cons  stantly  smaller  or  larger  than  the 
recorded  hydrographs,  the  error  may  be  in  the  rainfall  data  or  the  run¬ 
off  loss  parameters  used  in  HE C-l.  The  rainfall  data  should  be  exam¬ 
ined  and  possibly  the  rain  amounts  from  a  different  weighting  routine 
could  be  used  if  the  recorded  rainfall  showed  large  variations  between 
gages.  The  adjustment  in  loss  parameters  is  most  likely  more  effective 
in  changing  the  hydrograph  volume.  For  example,  changing  the  SCS  curve 
number  for  an  area  from  60  to  70,  the  runoff  from  a  lj-in.  rainfall 
would  increase  from  0.76  to  1.33  in.,  or  by  about  75  percent. 
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68.  Changing  the  streambed  or  overbank  roughness  factors  in  the 
HBC-2  input  will  affect  the  floodplain  storage  which  in  turn  will  modify 
the  hydrograph  attenuation  during  routing  by  HEC-1 .  Use  of  the  rough¬ 
ness  factors  for  hydrograph  adjustment,  however,  precludes  the  later  use 
of  roughness  factors  to  calibrate  the  stage-discharge  relationships  com¬ 
puted  by  HEC-2,  if  these  do  not  agree  with  the  relationship  between  re¬ 
corded  flood  peaks  and  observed  high  wutermarks. 

Choice  of  event  storms 

69.  The  expected  products  of  a  hydrologic-hydraulic  flood-level 
modeling  study  are  the  surface  water  profiles  and,  possibly,  the  flood 
damages  resulting  from  a  series  of  floods  of  specified  frequency  (re¬ 
currence  probability  or  return  period). 

70.  If  hydrologic  modeling  is  applied,  the  floods  must  be  gener¬ 
ated  from  appropriate  frequency  storms.  Even  though  it  is  recognized 
that  a  100-year  flood  event,  for  example,  may  be  produced  by  a  storm 
of  larger  or  smaller  return  period,  the  hydrologist  has  little  choice 
in  the  construction  of  his  model  but  to  assume  that  the  return  periods 
of  a  causative  storm  and  the  resulting  flood  are  identical.  The  basic 
rainfall  data  needed  as  hydrologic  input  are  then  the  magnitudes  of 
these  storms  of  specified  return  periods,  the  appropriate  rainfall  dura¬ 
tion,  and  time  intervals  to  be  employed  in  the  hydrograph  synthesis  and 
routing  processes. 

71.  The  set  of  return  periods  is  usually  specified  by  the  client 
or  the  agency  requesting  the  project.  Frequently,  however,  an  unnec¬ 
essarily  large  number  of  return  periods  is  specified,  resulting  in  un- 
managably  large  stacks  of  computer  output  which  do  more  to  confuse  than 
to  clarify  the  modeling  results.  The  requesting  agency  may  be  more  in¬ 
terested  in  only  examining  sumnary  tables  which  can  be  output  to  user’s 
requirements  or  may  choose  to  examine  alternative  methods  of  analyzing 
multiple  return  periods  based  on  data  obtained  from  a  prescribed  set  of 
return  periods.  For  example,  flood  peaks,  crest  elevations,  and  damages 
at  chosen  points  on  the  watershed  could  well  be  calculated  for  return 
periods  of  ?,  10,  25,  and  100  years  and  plotted  on  logarithmic  or  proba¬ 
bility  graph  papers  and  the  results  for  any  intermediate  return  period 
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read  from  the  graphs  with  as  much  reliability  as  when  the  modeling 
studies  are  repeated  for  all  return  periods.  Occasionally,  u  one  or 
even  one  half  year  flood  and/or  the  standard  project  flood  may  be  added 
to  the  above  list  of  four  basic  return  periods.  A  500-year  event  flood 
is  often  generated  for  flood  insurance  studies.  It  should  be  clearly 
understood,  though,  that  results  are  highly  suspect  because  rainfalls 
for  recurrence  intervals  exceeding  100  years  are  subject  to  large 
uncertainties.  It  is  the  responsibility  of  the  user  to  select  a  set  of 
return  periods  and  an  output  format  on  the  basis  of  cost,  time,  and 
level  of  accuracy. 

72.  The  rainfall  duration  chosen  should  be  at  least  as  long  as 
the  estimated  travel  time  through  the  watershed.  In  contrast  to  the 
Rational  formula  for  which  the  storm  duration  chosen  is  equal  to  the 
travel  time  and  applied  uniformly  not  only  over  the  drainage  area  but 
also  over  the  storm  duration,  the  rainfall  in  a  KEC-i  hyirograph  simu¬ 
lation  is  applied  over  up  to  1‘  0  time  interval:.,  usually  as  a  storm  of 
time-varying  intensity. 

73.  Storms  of  a  desired  return  period  and  duration  can  be  obtained 
from  S.  Weather  Bureau  Publication.  ’*  or  in  many  3tates  from 
more  detailed  maps  prepared  for  regional  rainfall  studies.  To  obtain 

a  time-varying  rainfall  input  it  is  suggested  that  hyetographs  be  con¬ 
structed  from  2«-hr  design  storms  using  'he  dimensionless  temporal 

12 

storm  distribution  developed  by  ICC ,  or  similar  rainfall  duration- 
intensity  curves  developed  by  other  agencies. 

7U.  If  hydrography  are  to  be  generated  for  largely  varying  areas, 

it  may  be  advisable  to  adjust  rainfall  as  a  function  of  watershed  area 

o 

as  prescribed  in  U.  S.  Weather  Bureau  TP  Uo  or  apply  an  elliptical 
storm  distribution  pattern  as  prescribed  by  the  U.  E.  Army  Engineers1" 
for  the  standard  project  storm.  If  the  inverse  distance  rainfall  weight¬ 
ing  program  (Appendix  C)  is  used  in  the  model  calibration  with  recorded 
hyetographs,  the  elliptical  storm  distribution  can  be  superimposed  on 
the  watershed  map  and  a  rainfall  assigned  to  each  rain  gage,  which  is 
then  converted  to  a  hyetograph  for  each  gage,  llie  weighting  program 
can  then  be  used  to  compute  a  weighted  hyetograph  for  each  minibasin. 


Generation  of  event  hydrograph 

75.  After  adjustment  of  watershed  parameters  by  model  calibration 
and  the  construction  of  input  hyetographs  by  procedures  described  ear¬ 
lier  under  "Choice  of  Eveut  Storms"  associated  event  hydrographs  can  be 
generated  with  the  use  of  the  HEC-1  program.  As  described  earlier,  the 
storage-outflow  relationships  for  the  Y-card  input  may  have  to  be 
generated  by  a  preliminary  run  of  HEC-2,  but  this  will  probably  be  done 
already  during  the  model  calibration  process. 

76.  The  HEC-1  input  deck  supplied  by  HYPPAR2  and  HEC-2  should  be 
checked  and  supplemented  where  necessary,  using  the  most  up-to-date 
HEC-1  input  instructions. 

77.  Upon  examining  the  HEC-1  output  hydrograph,  the  earlier  sug¬ 
gestion  to  keep  the  number  of  modeled  return  p>eriods  small  may  begin 

to  make  sense.  The  HEC-1  output  is  rather  voluminous  even  for  Just  or.e 
subbasin  with  1*0  minibasins;  the  output  for  a  large  watershed  and  10 
return  periods  could  fill  an  entire  room.  The  HEC-1  option  to  print 
only  the  most  essential  information,  namely  the  flood  peaks,  may  be 
useful  in  reducing  output  volume  as  well  as  providing  the  input  to  the 
final  HEC-2  run. 

Backwater  computations 

78.  Backwater,  or  flood  profile  computations,  car.  be  performed 
by  the  HEC-2  or  a  similar  program.  The  program  XSEC-PROC  can  provide 
the  entire  input  deck  except  for  the  title  cards,  the  J  cards  for  ini¬ 
tiating  the  program,  and  certain  options,  and  the  Job  ending  cards  Ed 
and  ER.  Some  special  optional  cards  such  an  X1*  or  C  cards  will  not  be 
provided  by  XSEC-PROC  either. 

79.  A  fair  amount  cf  visual  checking  and  comparing  of  the  HEC-2 
output  is  necessary  to  remove  occasional  inconsistencies  like  decreas¬ 
ing  flood  levels  with  increasing  discharges.  The  most  up-to-date  HEC-2 
version  and  input  instructions  should  be  used. 

Areal  delineation  of  floods 

80.  The  flood  levels  computed  by  HEC-2  must  be  extended  laterally 
to  plot  the  flood  boundaries.  As  long  as  the  stream,  as  well  as  the 
topographic  contour  lines  are  reasonably  ssieoth,  the  flood  level 


3U 


boundary  mapping  3houId  be  a  relatively  easy  task,  even  for  a  computer¬ 
ized  plotting  process  used  in  conjunction  with  a  digitized  topographic 
data  base.  Considerable  technical  Judgment  is  necessary,  however,  to 
plot  flood  level  boundaries  along  channel  bends  and  across  road  or 
railroad  embankments. 

8l.  For  two  separate  reasons,  it  is  strongly  suggested  that  the 
energy  line  be  used  instead  of  the  computed  water  surface  for  flood 
plotting.  First,  the  computed  water  surface  elevation  fluctuates  with 
velocity,  tending  to  be  low  in  swift  flow  reaches  of  the  stream  and 
high  in  slow  flow  reaches.  These  undulations  do  not  properly  represent 
overbank  conditions  where  water  surface  profiles  tend  to  be  much 
smoother.  The  computed  energy  line  changes  much  more  gradually  than 
the  computed  water  surface.  Second,  the  computed  water  surface  repre¬ 
sents  the  average  water  surface  in  the  stream,  disregarding  any  waves 
which  may  be  created  by  the  stream  velocity  and  carried  over  into  the 
floodplain.  The  energy  line  represents  the  average  water  surface  <  •  - 
vations  plus  the  kinetic  energy  or  velocity  head.  The  energy  line  thus 
constitutes  the  upper  limit  of  the  wave  height  which  can  be  generated 
ir.  the  stream,  and  sent  out  across  the  overbank3. 

Flood  damage  computations 

02.  The  estimate  of  potential  flood  damages  is  cften  the  most 
difficult  and  controversial  task  in  a  flood  study.  Damages  suffered 
during  previous  floods  may  depend  more  on  the  amount  of  advance  warning 
than  on  the  stage  of  the  flood  crest.  It  is  also  very  difficult  to 
estimate  the  relationship  between  flood  depth  and  damages,  particularly 
in  commercial  or  industrial  buildings.  Managers  or  owners  of 
industrial-commercial  enterprises  are  sometimes  reluctant  to  talk  about 
flood  damages  suffered,  or  they  may  exaggerate  the  estimates  in  order 
to  improve  the  likelihood  of  having  flood-control  projects  approved. 

03.  If  several  substantial  floods  have  been  experienced  over  the 
pant  few  years,  a  flood  discharge- flood  damage  curve  can  be  generated. 
In  a  systematic  modeling  project,  discharge  damage  would  best  be  tabu¬ 
lated  or  drawn  for  each  minibasin  separately. 

8h.  Frequently,  however,  too  few  floods  have  been  experienced  to 
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establish  a  reliable  curve,  and  potential  damages  are  estimated  by  tab¬ 
ulating  the  cumulative  value  of  perishable  property  located  at  or  below 
given  topographic  elevations.  This  method  results  in  a  curve  of  cumula¬ 
tive  damages  versus  flood  elevation  which  iB  not  very  useful  because 
the  peak  water  surface  elevation  of  any  flood  event  varies  with  its  lo¬ 
cation  relative  to  the  stream  (see  paragraph  8l).  Thus  the  potential 
flood  estimates  should  be  computed  individually  by  grids  small  enough 
that  the  flood  level  within  the  grid  may  be  considered  horizontal. 
Damages  can  then  be  related  to  discharges  by  interpolating  between  water 
surface  elevations  or  energy  grade  lines  at  the  nearest  upstream  and 
downstream  cross  sections  computed  with  HEC-2.  The  damages  resulting 
from  a  given  flood  may  finally  be  added  for  all  blocks  within  a 
minibasin. 

85.  In  conjunction  with  flood  frequency  analyses  the  discharge- 
damage  relationships  can  be  converted  into  curves  of  damage  versus  flood 
exceedence  probability,  as  shown  in  Figure  13.  Total  annual  expected 
damages  are  determined  by  integrating  the  damages  of  the  full  exceed¬ 
ence  probability  range.  Figure  13  illustrates  this  integration  graphi¬ 
cally.  The  total  damages  TO 

P 

TO  •  D  dP  (5) 

are  represented  by  the  area  under  the  damage-probability  curve.  In 
Equation  5,  D  i?  the  damage  caused  by  a  flood  of  exceedence  probabil¬ 
ity  P  and  P  is  the  exceedence  probability  of  the  largest  flood 
which  causes  zero  damages. 

86.  In  Figure  13  it  will  be  noticed  that  the  flood  damages  cor¬ 
responding  to  an  exceedence  probability  close  to  zero  are  large  and  un¬ 
known.  This  is  a  cosmon  dilemma  in  flood  damage  analyses.  Flood  peaks 
for  500-  or  1000-year  return  periods  could  be  computed,  but  results 
would  be  highly  uncertain;  a  manual  curve  extension  could  be  made,  but 
it  would  probably  be  equally  unreliable  (see  paragraph  71 ) .  Thus,  both 


exceedence  probability  and  damage  estimates  are  highly  unreliable  for 
extremely  large  floods.  This,  in  turn,  also  renders  the  determination 
of  average  annual  flood  damages  a  highly  speculative  process. 
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Excaadenc*  Probability,  Parcant 

Figure  13.  Damage-probability  curve  generated  from  f  1  ood-danage 
and  flood-probabi llty  curves 
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APPENDIX  A:  COMPUTER  PROGRAM  XSEC-PROC 


1.  A  computer  program  vats  written  to  convert  data  from  a  cross- 
section  data  file  to  standard  KEC-2  input  cards.  The  program  as  written 
takes  its  input  from  a  tape  and  produces  punched  cards  as  output. 

Simple  control  instructions,  however,  can  be  UBed  to  change  the  input 
and/or  output  type.  The  processes  followed  in  the  XSEC-PROC  program 
are  s^ovn  in  Figure  Al. 

2.  The  advantages  of  using  the  computer  program  XSEC-PROC  over 
conventional  procedures  are: 

a.  Field  cross-section  data  (i.e.,  rod  readings)  can  easily 
be  converted  to  HEC-2  input  data  deckB  (i.e.,  BT  and  GR 
cards)  provided  field  crews  are  avar-  of  some  simple  coding 
criteria. 

b.  Coding  for  either  normal  or  special  bridge  routines  is 
done  by  the  program;  the  operator  has  only  to  decide  which 
routine  is  most  applicable  for  a  given  situation. 

c.  Many  of  the  variables  required  on  the  SB  card  (i.e., 
trapezoidal  area,  bottom  width,  side  slaps,  etc.)  used 

in  the  special  bridge  routine  are  computed  within  the  pro¬ 
gram  and  reduce  the  time  required  for  these  factors  to  be 
determined  manually. 

d .  Repeated  sections  ewe  automatically  raised  or  lowered  de¬ 
pending  of  the  slope  of  the  stream  bottom  between  two 
measured  sections. 

£.  The  probability  of  human  error  in  coding  HEC-2  data  deckB 
is  greatly  reduced. 

3.  The  input  to  XSEC-PROC  consists  of  the  regular  cross-section 
survey  data,  preceded  by  supplemental  data  lines  labeled  A  to  K  as 
described  below.  The  identification  letters,  entered  in  column  1  of 
each  line,  are  not  mandatory,  but  they  make  it  easier  to  read  the  coded 
information. 


Line  A.  -  Title  or  Comsent  Line.  Format  (l7Ak,  IU) 


Columns 


69-72 


Variable  Value  _ Description _ 

COMM  (I)  Any  alphanumeric  information. 

IC  0  or  1  If  1C»1,  read  another  A  line. 


Al 


PHASE  I  INPUT  DATA  RE  AOING 


READ  NUN  TIM  I  AND  COMMENTS  ON 
A  MNES 


RE  AO  SECTION  NUMBER.  LOCATION 
NUMBER  Of  SECTION  TRANSEERS  ANO 
INOE*  EOR  BRIDGE  IBR  ON  0  TINE 


NS  TO  -  NO  Of  OOANSTREAM  TRANSEERS 
NSTU  NO  Of  UPSTREAM  TRANSEERS 
If  IBR  0  NSEC  *  NS  TO  ♦  NSTU  •  I 
If  IBR  0  NSFC  *  NSTO  •  NSTU  •  1 


REPEAT  NSfC  EOOPS 


RE  AO  TRANSFER  DISTANCES.  NUMBER  Of  Nf* 
N  ANO  OR  Q  VALUES  NN  AND  NO  AND 

INOE*  1*3  fOR  LEVEES  ON  f  LINE 


READ  NC  CARO 
ENTRIES  ON  f  UNI 


RE  AO  NH  CARO 
ENTRIES  ON  G  LINE 


TO 

SHEET  ) 


Figure  Al.  Flov  chart  of  XPEC-PfTOC  (sheet  1  of  M 


MM 


at  AO  XI  CARO 
(NT  HIES  ON  I  lINf 


HAVE  NSEC  SETS  Of 
f  TO  I  UNIS 
BEEN  Rf  AO’ 


Ht  AO  BRIDGE  DATA 
J.X  HNE5 


RE  AO  COOEO  SURVEY  DATA 
ST  A  Ht  ROO  READING  IROGI 
ANO  COOES  AA 


RMASE  1  OAT  A  PROCESSING 


Figure  A 1  (sheet  2  of  M 


PHASE?  OW  A  PROCESSING 


SEPAR  AT f  SURVEY  DATA  INTO 


BRIDGE  SECTION  DATA 


GROUND  SECTION  DATA 


CHECK  TOR  INCREASING  STATION  SEQUENCE 
ANO  CORRECT  If  NE CESSARY 


SEARCH  EOR  BEGINNING  ANO  ENO 
Of  BRIOGE  OPENINGS 


PREPARE  ARRAT  Of 

GROUNO  SECTION  STATIONS 
ANO  ELEVATIONS 


SPECIAL  BRIOGE  ROUTINE 
REGISTER  BT  CARO  OATA  WITHIN 
SPECIFIED  WEIR  LIMIT  STATIONS, 
UNLESS  LENGTH  Of  HORIZONTAL 
WEIR  IS  PREDETERMINE 0 
COMPUTE  BRIOGE  OPENING  AREA 
AND  TOTAL  PIER  WIDTH 


NORMAL  BRIOGE  ROUTINE 
REGISTER  BT  CARO  DATA  WITHIN 
BRIOGE  OPENING  LIMITS  WHEREVER 
LOW  CHORD  ANO  TOP  Of  ROADWAY 
SURVEY  DATA  POINTS  WERE  TAKEN 
MODIFY  STREAM  BEO  SECTION  BY 
INCLUDING  PIERS  AS  SOHO  BOUNDARIES 


COMPUTE  REACH  LENGTH  ANO 
STREAM  8EO  SLOPE  BE  TWEE N  SECTIONS 


PUNCH  (OR  TAPE  I  Hf  C  ?  INPUT  STARTING  WITH  LOWEST 
TRANSFER  SECTION  WHERE  APPLICABLE 


PUNCH  QT  CARO  If  NO  >  0 
NC  CARO  If  NN  ■  3 
NH  CARD  If  NN  >3 


Fif^ure  A1  (sheet  !*  of  l») 
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Line  B.  -  Printing  and  Units  Codes. 


Columns 


1 

2-h 


5-8 


Variable 


A  (2) 
IPRI 


IUN 


Value 


B 

0  or  ♦ 


0,  1,  or  2 


Format  (AI,  13,  lL) 
_ Description 


Printing  instruction.  If  1PRI>0, 
XSEC-PROC  input  data  are  printed 
out. 


Measurement  units  index  for  dis¬ 
tances  between  sections  and  water 
shed  outlet. 


If  IUN=0  distance  is  in  feet. 
If  IU?«=1  distance  is  in  miles. 


If  IUNB2  distance  is  in  kilometres 
and  all  other  dimensions  in 
metres. 


Lines  C  to  K,  where  applicable,  and  survey  data  lines  are  repeated 
for  each  section. 

Line  C.  -  Initializing  and  end-of-file  line.  Format  (Al,  812,  Lx,  A3, 


2LF2.0) 


Columns 

Variable 

Description 

1 

A  (3) 

C 

2-21 

Mi  sc. 

blank 

Various  integers  which  must  be  set 
to  zero  prior  to  processing  a  new 
section. 

22-2L 

AA 

blank 

Read  section  data. 

EOF 

End  of  file  is  reached.  Exit. 

25-72 

Mi  sc. 

blank 

Various  real  variables  which  must 
be  set  to  zero  prior  to  process¬ 
ing  a  new  section. 

Line  D. 

-  Teneral  Section 

Information.  Format  (Al,  17,  LlL,  LF6.0,  23X, 

ID 

Columns 

Variable 

Value 

De script  ion 

1 

A  (J») 

D 

2-8 

ISEC  (J) 

♦ 

Cross-section  identification  number. 

Identification  number  of  first 
cross  section  of  a  tributary. 
ISEC  must  be  the  same  as  I SEC  of 
cross  section  Just  below  the 
Junction. 


Line  D. 

(Continued) 

Columns 

Variable 

Value 

Description 

9-12 

NSTD 

♦ 

Number  of  copied  cross  sections 
transfered  downstream.  Minimum 
of  1  at  a  bridge. 

13-16 

SSTU 

♦ 

Number  of  copied  cross  sections 
transfered  upstream.  Minimum  of 

2  at  a  bridge. 

17-20 

IBR 

0 

No  bridge  at  this  section. 

1 

Normal  bridge. 

2 

Special  bridge. 

21-2L 

ITR 

0 

1 

At  all  times,  except  if  last  trans¬ 
fered  cross  section  is  located 
Just  below  a  tributary.  ISEC 
will  be  increased  for  this  trans¬ 
fer  section. 

25-30 

XRIV 

♦ 

Distance  from  section  to  watershed 
outlet,  in  feet,  miles,  or  kilo¬ 
metres,  depending  on  the  value  of 
I UN  on  line  B. 

31-36 

37-1*2 

STCH  (1) 

STCH  (2) 

♦ 

Left  and  right  channel  boundary 
stations. 

l»3-!»8 

THETA 

Angle  by  which  the  section  deviates 
from  the  normal  to  the  channel , 
in  degrees. 

72 

IC 

0  or  1 

If  IOi  an  A-Line  is  read  after 
this  line,  containing  verbal 
descriptions. 

Lines,  E,  F,  G 

,  H,  and  I  are 

repeated,  where  applicable,  a  number 

of  times 

equal  NSTD 

♦  NSTO  ♦  1. 

Line  E. 

Transfer  Section  Line. 

Format  (Al,  13,  21 A,  5F6.0) 

Columns 

Variable 

Value 

Description 

1 

A  (5) 

E 

2-U 

NN  (I) 

0 

No  new  n  values  for  transfer 

section  I. 


3  3  new  n  values  and  possibly  new 

contraction  and  expansion  coef¬ 
ficients  will  be  read  on  F  line 
for  NC  card. 


A7 


Line  E.  (Continued) 

_ Description _ 

NN  new  values  will  be  read  on  G 
line  for  NH  card. 

No  QT  card  needed. 

NQ  discharge  values  to  be  read  on 
H  line  for  QT  card. 

No  X3  card  needed. 

Standard  X3  card  to  be  punched. 

Levee  values  will  be  read  on  I  line 
for  X3  card. 

Distance  of  transfer  section  from 
surveyed  section  (-  for  down¬ 
stream,  ♦  for  upstream  transfer 
section) . 

Drop  or  lift  of  transfer  section  I 
with  respect  to  surveyed  section, 
for  entry  in  field  9  of  XI  card. 
If  DYTR  «  0,  drop  or  lift  will  be 
computed  fro®  slope  and  reach 
length.  DYTR  must  not  be  0  in 
first  section. 

No  section  width  multiplier  needed. 

Section  width  multiplier  for  field 
8  of  XI  card. 

Overbank  reach  length  equals 
channel  length. 

Left  ar.d  right  overbank  reach 
length  increments,  respectively. 

Line  F.  (If  NN  (I)  *  3)  -  NC  Card  Entries.  Format  (A 1,  F5.0,  11F6.0) 


Columns 

Variable 

Value 

Description 

1 

A  (6) 

F 

2-6 

VNC  (I,  II) 

0  or  ♦  Entries 

for  NC  card. 

7-12 

II  -  1,  5 

13-18 

etc. 

Line  G. 

(If  NN  (I)  >  3) 

-  NH  Card  Entries. 

Format  (Al,  F5.0,  11F6.0) 

Columns 

Variable 

Value 

Description 

1 

A  (7) 

G 

Columns 

Variable 

Value 

>3 

5-6 

NQ  (I) 

0 

♦ 

7-8 

1X3  (I) 

0 

1 

2 

9-lU 

DXTR  (I) 

-  or  * 

15-20 

DYTR  (I) 

-  or  ♦ 

21-26 

XV 

(I) 

0 

♦ 

27-32 

DXL 

(I) 

0 

33-38 

DXR 

( I ) 

♦  or  - 

A8 


Line  G. 

(Continued) 

Columns 

Variable 

Value 

Description 

2-6 

VNH  (I,  1) 

4 

First  n  value  for  NH  card. 

7-12 

XNH  (I,  1) 

4 

End  station  of  first  n  value,  for 

NH  card. 

13-18 

vnh  (i,  2) 

4 

Second  n  value. 

19-2U 

XNH  (I,  2) 

4 

End  station  of  second  n  value. 

etc. 

etc. 

* 

etc. 

repeat  NN  { i ) 

times 

Line  H. 

(If  NQ  (I)  >  0)  -  QT  Card  Entries.  Format  (Al,  F5.0,  11F6.0) 

Columns 

Variable 

Value 

Description 

1 

A  (8) 

2-6 

QT  (I,  1) 

4 

Flow  1  for  QT  card  in  Auxiliary 
Section  I. 

7-12 

QT  (I,  2) 

Flow  2  for  QT  card  in  Auxiliary 
Section  I. 

13-18 

QT  (I,  3) 

Flow  3  for  QT  card  in  Auxiliary 
Section  I. 

etc. 

etc. 

etc. 

Line  I. 

(If  1X3  *  2)  - 

X3  Card  Entries.  Format  (Al,  F7.0,  8F8.0) 

Columns 

Variable 

Value 

Description 

1 

A  (9) 

I 

2-8 

9-16 

17-2b 

X3  (I,  1) 

X3  (I,  2) 

X3  (I,  3) 

0  or  ♦ 

KEC-2  Entries  in  X3  cards  for  levee 
and  encroachment  data. 

etc. 

etc. 

Line  J. 

(If  IBB  >  0)  - 

Bridge  Data.  Format  (Al,  F5.0,  11F6.0) 

Columns 

Variable 

Value 

Description 

1 

A  (10) 

J 

2-6 

BRV 

4 

Bridge  width,  in  feet  or  metres. 

7-12 

PV 

4 

Pier  width,  in  feet  or  metres. 

13-18 

SHR 

4 

Solid  height  of  railings. 

19-2>» 

STCHB  (1) 

4 

Left  channel  bank  station  in  bridge 

section. 
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Line  J.  (Continued) 


Columns 

Variable 

Value 

Description 

25-30 

STCHB  (2) 

■f 

Right  channel  bank  station  in 
bridge  section. 

Line  K. 

(If  IBR  « 

2)  - 

Additional  Special  Bridge  Data.  Format  (Al, 

F5.0,  11F6.0) 

Columns 

Variable 

Value 

Description 

1 

A  (11) 

K 

2-6 

7-12 

13-18 

SB  (1) 

SB  (2) 

SB  (3) 

♦ 

3  special  bridge  coefficients  for 
fields  1-3  in  SB  card.  (See 

HEC-2  manual  for  suggested  values. 

19-21* 

WL 

0 

Weir  length  will  be  determined  by 
limits  WSTL,  WSTR  of  BT  cards. 

♦ 

Special  bridge  effective  weir 
length.  Weir  should  be 
horizontal . 

25-30 

SS 

0  or  ♦ 

Abutment  side  slope  jr  . 

31-36 

WSTL* 

♦ 

Starting  station  of  weir  described 
by  BT  cards. 

37-42 

WSTR* 

■f 

Ending  station  of  weir  described 
by  BT  cards. 

•  WSTL 

and  WSTR  should 

be  left  blank  if  WL  >  0,  and  vice  versa. 

Cross-Section  Survey  Data  -  Format 

(3F6.0,  3X,  A3,  2F6.0) 

These  lines  do  not 

hove 

identification  letter  in  column  1. 

Columns 

Variable 

Value 

Description 

1-6 

ST 

0  or  * 

Station,  in  feet  or  metres,  of  a 
survey  observation. 

7-12 

HJ 

0 

Previous  height  of  instrument  line 
of  sight  will  be  used. 

♦ 

Sew  height  of  instrument  line  of 
sight. 

13-18 

RDG 

♦ 

Rod  reading  at  ground  or  road. 

22-24 

AA* 

A 

Alphabetic  code  word,  to  be  matched 
against  registered  codes  which 
direct  the  program  into  different 
branches  of  operations.  See 
below  for  code  words. 

Cross-Section  Survey  Data  (Continued) 


Columns 

Variable 

Value 

Description 

25-30 

YLCD 

■f 

Rod  reading  at  low  chord  of  bridge. 

31-36 

KURD 

♦ 

Rod  reading  at  top  of  roadway 
above  bridge  opening. 

*  Codes  CODE  (i)  against  which  readings  AA  will  be  matched. 


i_  CODE  (i) 

1  3  blanks 

2  BPR 

3  BSB 

l  CLR 

5  EBR 

■f 

'6.  END 

7  EWS 

i 

8  EOF 


_ Description 

Normal  ground  station. 

Bridge  pier  station. 

Streambed  station  under  bridge,  no  pier. 

Center  line  of  road  outside  bridge. 

Edge  cf  bridge.  More  than  one  bridge  opening  is  accept¬ 
able,  and  an  EBR  line  should  be  coded  at  the  beginning 
and  end  of  each  bridge. 

End  of  section  stations. 

Edge  of  water  surface  at  normal  (fairly  dry -weather  flow). 
Only  one  edge  of  the  stream  needs  to  be  coded  and  the 
reading  will  be  used  to  estimate  normal  stream  slope. 

End  of  file.  Causes  program  to  terminate  data  reading 
and  processing. 


Definition  of  Variables 


A  (I) 

AA 

Al,  A  (2), 
A  (3) 

ABO 

B1 

BT  (I,  II) 
BT2,  BT3 


*  Array  of  identification  letters  (from  A  to  J)  to  be 
entered  into  Column  1  of  XSEC-2  input. 

=  Alphabetic  variable  to  be  read  as  code  from  cross- 
section  data  lines  in  input,  columns  22  to  2h. 

*  Arrays  of  1,  2,  and  3  constants  set  to  zero,  to  be 
called  on  when  wanted  for  output. 

*  Area  of  bridge  opening,  computed  internally  from  the 
survey  data  by  the  trapezoidal  rule. 

=  Constant  set  equal  to  1,  to  be  called  on  when  wanted 
for  output. 

*  Matrix  of  output  values  for  BT  cards,  for  up  to  50 
stations  per  section. 

=  Variables  for  temporary  assignment  of  BT  (i,  2)  and 
BT  (1,  3)  values. 


All 


CODE  (I) 


COMM  (I) 

DXL  (I) 
DXR  (I) 
DXTR  (I) 

DYTR  (I) 

DY,  DZ 


EL  (I) 
ELB  (I) 
HI 
HJ 

IBR 

IBT 

I  SB 


*  Alphabetic  code  entered  in  data  block.  Variables  AA 
read  from  tape  are  compared  with  CODE  (i)  to  identify 
cross-section  stations  as  roadbed,  bridge  pier,  and 
other  features.  CODE  (i)  from  i  »  10  to  22  is  used  for 
writing  HEC-2  codes  XI,  GR,  etc.,  and  for  printing  out 
messages . 

*  Array  of  17  comment  words  (alphabetic)  read  in  subrou¬ 
tine  COUNT  and  printed  out  but  not  punched  into  HEC-2 
input . 

Array  of  distances  by  which  channel  reach  lengths  should 
be  lengthened  or  shortened  for  left  and  right  overbanks. 

*  Array’  of  distances  (-  =  downstream,  ♦  *  upstream)  by 
which  surveyed  cross  6ection6  should  be  transferred  to 
create  artificial  sections. 

*  Array  of  y-increments  by  wnich  a  transferred  cross 
section  should  be  moved  up  or  down  relative  to  surveyed 
cross  section. 

*  Left  and  right  height  of  each  trapeioidal  section  of 
bridge  opening.  Used  to  compute  area  ABO  of  bridge 


STB  STE 

Figure  A2.  Bridge  opening  representation 


*  Array  of  elevations  for  GR  cards  of  ground  sections. 

®  Array  of  elevations  for  GR  cards  of  bridge  sections. 

*  Surveying  height  of  instrument. 

*  HI  recorded  in  surveying  notes.  As  long  as  HJ  *  0, 
previous  HI  prevails.  When  HJ  >  0,  HI  is  set  equal 
to  HJ . 

*  Bridge  type  identification  (l  ®  normal,  2  *  special)  in 
input  data. 

*  Counter  to  indicate  whether  a  station  lies  within  limits 
of  BT  statements  in  special  bridge  routine. 

«  Counter  for  bridge  abutment  stations. 
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10,  110 


I SEC  (I) 


1X3  (I) 


NN  (I) 


HQ  (I) 


Constants  set  equal  to  0  and  10,  respectively,  for  use 
in  output. 

Do  Loof  parameter  for  total  number  of  transfer  sections, 
and  several  other  do  loops. 

Minor  Do  Loop  parameter. 

Indication  that  input  statement  will  be  followed  by  a 
comment  statement  (if  IC  >  0). 

I.D.  number  of  cross-section  I. 

Indicator  to  send  routine  bacfc  to  statement  to  gener¬ 
ate  new  NC  or  NH  cards  before  or  after  bridge  opening. 

Counter  for  stations  inside  bridge  openings. 

Instruction  to  print  back  inpqt  data  if  IPRI  >  0. 

Temporary  section  I.D.  number.  If  an  ISEC  <  0  is  found, 
indicating  the  start  of  a  tributary,  ISCN  is  set  ■  -ISEC 
so  that  the  corresponding  downstream  section  of  I.D. 

No.  ISCN  car.  be  located. 

Indicator  that  a  tributary  Joins  above  an  upstream  trans¬ 
fer  of  this  section. 

Entry  for  field  1  in  XI  card. 

Instruction  to  produce  an  curd  for  transfer  Bection  I. 
If  1X3  ■  1  produce  standard  X3  card.  If  1X3  *  2  produce 
X3  card  as  read  from  input  tape. 

Instruction  for  distance  conversion  from  miles  to  feet 
if  I UN  =  1,  kilometres  to  metres  if  IUN  *  2. 

Counter  for  cross  sections.  Used  as  subscript  for  sec¬ 
tion  data. 

Counter  for  stations  in  ground  section. 

Counter  for  stations  in  bridge  sections. 

Counter  for  BT  card  stations. 

Number  of  new  roughness  n  values  (if  any)  for  transfer 
section  I  to  be  punched  onto  NC  or  NH  cards.  New  n 
values  may  be  dictated  for  up  to  10  transfer  sections 
per  surveyed  section. 

Temporary  substitute  for  NN  (I). 

Number  of  Q  values  to  be  read  at  transfer  section  I  from 
data  and  punched  onto  QT  card. 

Temporary  substitute  for  NQ  (I). 

Number  of  downstream  transfer  sections. 

Number  of  upstream  transfer  sections. 
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NSEC 

PV 

QT  (I, II) 

RDG 
RL  (I) 

RVRL 

Rl,  R2,  R3 
SB  (I) 

SHR 
SL  (J) 

SS 

ST 

STA  (X) 

STAB  OCX) 
STB,  STE 

STCH  (1), 
STCH  (2) 
STCHB  (I), 
STCHB  (2) 
THETA 

VNC  (I,  II), 
VWH  (I,  II) 

VBO 

WL 

VPR 

WS 

WSTL,  WSTR 

XEB 

XKB 

XRN 


NSTD  ♦  NSTU  ♦  1  »  total  number  of  XI  cards  to  be  gener 
ated  for  cross-section  J. 

Pier  width. 

Flows  to  be  entered  into  QT  card . 

Rod  reading  in  survey  notes. 

Reach  lengths  for  left  bank,  right  bank,  and  channel. 
Bridge  deck  rod  reading. 

Standard  reach  lengths  used  on  bridge  sections. 

Entries  for  SB  card. 

Solid  height  of  bridge  railing. 

Stream  slope  at  section  J. 

Side  slopes  of  bridge  abutments  (for  SB  card). 

Station  distances  in  feet  (or  metres)  read  from  survey 
data. 

ST  value  at  station  K  in  ground  section. 

S.  value  at  station  KX  in  ground  section. 

Beginning  and  ending  stations  of  trapezoidal  sector  or 
bridge  opening.  See  Figure  A2. 

Left  and  right  channel  stations,  for  XI  card. 

left  and  right  channel  stations  for  bridge  section. 

Angle  of  section  and  with  respect  to  a  normal  stream. 
Entries  In  field  II  of  NC  or  NR  cards. 

Width  of  bridge  opening  (entry  for  SB  (5)). 

Weir  length  (entry  for  SB  (M). 

Total  pier  width  (entry  for  SB  (6)). 

Water  surface  elevation  at  normal  stream  flow. 

Left  and  right  stations  to  be  used  for  weir  length  in 
special  bridge  routing  if  WL  ■  0  and  8T  cards  are  used 

Indicator  of  whether  a  station  is  inside  or  outside  a 
bridge  opening. 

Floating  point  definition  of  KBT. 

Floating  point  definition  of  NTi. 


XQT 

XNH  (I,  J) 
XRIV 
XR  (J) 

XXR  (I) 
XXS,  XX9 

X3  (I,  II) 
X8,  X9 
YY8,  YY9 
Y  (IOP) 
YLCD 


■  Floating  point  definition  of  HQT. 

■  Station  in  NH  card. 

=  River  mile  or  Kilometres  from  survey  notes  for  section  J. 
=  XRIV  converted  to  feet  or  metres. 
a  XR  values  for  transfer  sections. 

*  Entries  for  fields  8  and  9  in  X3  card  downstream  from 
bridge . 

*  Entries  for  X3  cards,  if  specified  in  input  data. 

*  Entries  for  fields  8  and  9  in  XI  card. 

*  Same  as  XX 8 ,  XX9,  but  upstream  from  bridge. 

■  Low  chord  elevation  in  bridge  .opening. 

■  Low  chord  rod  reading. 
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XSBC-PROC  Listing. 


c  MKm>XAM  KSLC-MMUf  II)  iNfMJl  lll-CK  Mix  M(-C -? 

DlMCNMUN  A  I  1  ?  I  ,  A3I3I.  MHnO.31%  U«J*  <  Mil  .  HXUlU).  OXHllO). 
l  nkiMiioi.  nnKiim.  t-iiiuoi,  llmioid,  i^nvmi.  lAtuoi,  nniioi 
P  «  NO  I  1 1 )  I  «  OlllO.lSI,  KLIi).  «1I3).  K?<  ft  ,  Nil)).  SMI  l  01  .  SU  *>1101  , 

J  SIAI100).  SI  AM  ||)0)  ,  SICM(?I,  NTCHMIP),  VNUlfl.bl.  VNHllD.POI. 

-  XNM  I  1 1)  •  ?()  )  •  XMISOO),  *010),  UKIIOIII,  K?<  |0I  .  X3II0.V),  V  I  JODI  , 

*>  f«Slt>00) 

OAIA  CUOt-  ,  3HHHK  ,  3MM  SK  .  )M<  L*<  •  3Hf- MK  .  3MLND  .  3ML  WS  .  )Mt  <  D- ,  3MI1UM. 

I  PMNC  .  I  ,  <>M*  1  ,  ?HK?  ,  ?MX  3,  PUNK  ,  ?MH  I  .  ?MSN  ,  AMt.KII,  jMUNl  .  iMhK  |  ,  iMI) 

PL»t  ,  IM  M.AMAfcl.lH  Ml,  Antes,  AM  Kl.JMUIM.AM  ,  AM  / 
c  Ml  Al)  Mx  I N I  .MUNCH  I W  JAM  INStKUCt.  AND  lINlI  IIMlUNS 

C  IDN»0  -  ►M-l.  I  HN*  1  -  MILLS.  |UN*P  -  KJLIIMtltKS 
MLAD  I  ?  •  VO  I  A  I  1  )  ,  IHKI.IDN 
I  M  I  KK  |  .f.I  .0  )  MX  IN  I  wv 
H  •  ?>  ♦  ?*|UN 
1/  <  /A  ♦  ?  ° | UN 

1M  I  MK  I  .01  .«)  I  MM  I  Nl  MV,  All),  IfK  I  .  I  UN, CllUt  ( 1. 1  I  •  (  DDL  1 1.  P  I 
MV  MIKMlII  •  *  ,A|,?tv,<  CKIlSS  SfcCIlDN  LOCAIl'IN  HtASHKH)  IN'.PA)) 

CALL  CUMMNI ( I I | 

J  »  0 
10  «  0 
I  10  ■  |0 
?0  J  »  .l*t 

C  M  t  AD  ML  ana  VAX  1  AMLtS 

HtAD  t?,v|l  At  3  1  ,X  tKK.KHf  <  IH  I  ,  It  M.  AA.ANO.  Mt»H  .ML  «  tfrtl),  A3 
I f I AA. tD. CUOt l« I IGU  ID  t>l 
Dll  l  I  •  1,10 
*?<  I  )  •  0. 

DD  1  I  I ■ 1  .  V 
I  1311,111  •  0. 

C  ft  t  AD  StdlDN  NUMMtK  AND  CMAMACltHlSIlCS 

Ki-tii  (2.921  a  i  -  ) ,  I stCi.n  .  ns  Id,  NS  YU.  IHK,  I  YK.XHI  V.STCM.THflA.IC 
1M  |f*M  I  .t.r.O  |MM  Ini  101,  A  I  A  )  ,  IStCtJI.  nsid,  NSID,  Imk,  |  TM.  AM  I  V  , 
1 SICM. IMM  TA 

iHtTA  *  THt  1  A*  3.  lO  t>/ 1  MO. 

If- 1  IC.OI  ,0)CALL  CDmmnT  I  I  mm  I  I 
NSfcC  •  NSID  «  nSTd  ♦  1 
C  CMfcC*  f-UM  IKAMDfX  StCUUNS,  KtAU  DAIA 
mi  H  •  i ,Nstc 

Mf  AO  IP, Vi)  A  I  S  )  ,  NN II),  NO  I  I  I  , I  X  3  I  I  I ,  DXIMIII.  DVTMIII.  AMll), 
1I)*L  I  1  I  •  DAK  I  I  I 

IMI  f*K  I  .  1. 1  ,  ':  )  f*M  I  N I  l  OP, AIN)  •  NN  I  I  I  ,  Noll).  I  a  4  I  I  )  .  OX  IN  I  |  I  .  DVlMII 
I  I  .  *  Ml  |  I  ,  U*L  <11,  OXK  <  |  ) 

1  M  J.to.  DSL  I  I  )  »DY  IM  I  |  )  /UXlN  I  I  I 
IMXM  I  l.fcO.O.  I  A  N  I  |  1*1. 

NOT  •  NO  I  | I 
NN!  *  NNI 1  I 
J  f ( NNI -3 ) 17,IM,1V 

l  X  Nt  AD  IP. VO  A  I  A  I  .  I  VNC  (I,  111.  11*1. Ml 

1  f- I  IKK  |  ,f.t  .OIMM  INI  103,  AIM,  (  VNC  I  l  ,  1  I  1  ,  1  I  •  1  ,M 
Ul  TU  17 

IV  HLAD  I?. VO  A  I  7 )  ,  I  VNM (|,|l),  KNH 11,11),  11*1, NNtl 

TM  I**K  |  .  (.1 .0  IMKlNl  |  0  3  ,  A  I  /  I  ,  (  VNM  <|,l|l,  KNM II,  II).  11*1.  NN  I  I 
17  JMNOI  .MJ.OII.U  III  ?! 
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>  iv  .--i4 


Kt-Atl  <  ^  .  S>  ••  I  AIM.  (u  I  (  I  ,  |  II  *  1  1  *  l  .NO  I  I 

IM  It’**  I  .(.I  .OlPKlNl  l  •»<.  .  AIM,  (UT  II  .  Ill  .  I  I  ■  I  .NU!I 

7  1  M.n  III  ✓ 

Kt-An  I/.1M  AIVl,  (  A  3{  I  .  I  J  I  .1  I  *1  .V  \ 

IM  IPK  I  .(iT.lllPHlNI  I  IK.,  tl.l,  < X  3<  I  . I  I  I .  I  I  *  I .VI 
7  CUN  I  lNllt- 

Kt-AII  HK  |  !)(>►  1 1 A  I  A  II-  ANY 
I  t  (  |  KM  .Ml. 01  t.ll  III  40 

Kf-AO  17. VM  A  I  I  (1  I  .  HKH,  V  *  ,  SMK.  SltHH 
I  M  I  PH  |  .(.1  . O  I  HK  I  NI  UK..  MUU.HKH,  PH.  SMK  .  SICKH 
IM  IHK.MI.;|ilMI)  (7.VM  A  I  1  |  I .  (  SHI  I  I  I  .  |  I  ■  |  .  1 1  .  HI. ,  SS.HSTl. ,  uSI  K 

|H  IPH  |  .tO.  1  .  ANII.  IMK.MJ.7 |PK |Nl  1  04.  A  I  1  1  I  •  (  SH  I  I  I  •  I  ■  1 . 1 1  •  HI..SS 

IHSU.  WMH 

MM. IN  Kf  A 1 1  I  N(*  r.KIlSS  S  t-  C  T  I  I  IN  SlIHVf-V  IIA1A 
30  K  f-  Al  I  17. VM  SI.  M.l,  KIM.,  AA.  YLCO.  KHKII 

IM  |PK(  .(.I  .OlPKlNI  JOS.  SI,  MJ,  HIM..  AA  .  YI.CI).  KwHIl 
|  I-  (  M  l  .1.1  .1).  )M|  «M.|  , 

it-IAA.f-O.r.llUirlMli.ll  III  1 

I  ►  I  AA  ,t-  >  .emit  I  t  I  .!»*  .  AA.t-o.CUIIt  I  1  >11.11  III  40 
K1  «  K<*1 
S I  AM  IKK  I  ■  SI 

fiKiui  «  mi  -  mu, 
uu  rt  l  •  l .  J 
m  1 1  1 1  •  l . 

«? I  I  I  «  HKH 
"i  Kill)  »  HKH  ♦  7. 

HK  ANf.Mtlt- I-  MW  SMt-CIAL  HklOlilr 
I  M  I  MW  .(-0.7  lull  III  A«. 

IM  A  A  .  fcO.CUIIf-  I  J  I  .UK  .  AA.  tU.CUUM  A  I  jr.ll  III  Id 

I  M  KHKII.f-O.O.  H.II  Til  37 

KHl  •  KHl*l 

HllKHI.JI  »  SI 

hT(kmI,7I  ■  m|  -  mkhII  ♦  Smk 

HT IKHI .31  "Ml  -  V|.CO 

l»-(AA.Nr.cnoi-i7iii.ii  Tn  jo 

HI  (KHI  ,|l  •  SI  -  PH/7. 

37  ST  Aft  I  AA  I  •  ST  -  PH/7. 

KK  •  KK»1 

STahikki  «  sTahikk-ii 
H.HIKKI  «  M|  -  Vf.CH 
KK  a  KK  * | 

STAHIKKI*  STAHIKK-ll  »PH 
6UUS«I  «  CtHIKK-li 
KK  «  KK • 1 

SIAH(KK)  *  STAHIKK-ll 
f  C  H  I  KK  I  «  tf.HIKK-31 
M  TU  30 

34  IM  STAHIKK  I  .M'.HSIL  I  IH1  »1 
lMNi.(.i.n.  iuu  Tn  7H 
I  M  |  M  T  .t-O.OIf.ll  Til  30 
H  T  7  «  M.HIKKI 

HT  3  *  n. 

7«  I  m  a  a  .  mi.cuom  s  1 1 1  tn*  i  mi.  i 

X|-H  .  M.UAT  I  (t-Hl/7.  -  PL  IIAT  (  I  f-H/7  | 

1  M  M-H.t-O.O.  .  ANl).  A  A.  Nf  .  CUUk  I  M  I  1 1)*»  *  O 
IM  MrH.tO.O.ANll.AA.NF.CIinMM  IGIT  Til  3* 

I  H  A  A  .Ml.  CUlir  I  7  I  IHPK  *wpk  *ph 
KIP  *  |IIP  ♦  f 
STT-  *  ST 

VI  IIIPI  *  Ml  -  VLCll 


|  M  MWUH.1.1  .O.  )  VI«H1-KHKII*V<M 

I  M  UIH.Hi.n,  I  v  (  Iok  >  ay |  |n)>-  1 ) 

Y  M  a  Y  I  I  I  >►*  ) 

II?  *  VlllltM-tl,  H|KK) 

!►  10/  .1  I 

I  ►  »  Hip  ,M).  1  )i.n  III  )•> 

AMD  •  A M| |  ♦  (OYall?  )•<  S!t~StH  )/?. 

M  Ml  l  »  WMI)  ♦  Sit  -  SIM 
3">  AIM  ■  Alt 
IIV  »  11/ 

I  M  Ml.  .  (  I  .0.  )  Ml  III  30 
IMHMKII.tll.ll.ll.il  111  II 
Ml?  ■  Ml  -  HMHI)  ♦  SMK 
Mil  .  Ml  -  VI.  f.|) 

3s  I  ►  (  WI..I.I  .f)  .  I'M'  ID  30 

KMI  a  AMI  ♦  J 
MI  I  KM  I  . 1  )  *  A  I 
Mil KMI,?)  «  Ml? 

MI  I  KMT .  3 )  a  M  T  3 

IMI  ltM»  ||IM  I  .►!).  1  I XKM.HT (KMI , 3) 

I  *■  «  I  I  tH>*|l.“».MJ.  1  )YYM*M)(KMl  ,?) 

1 H  I  KMKI).(,T  .0. IXXV  a  MI (KMI , 31 
IMKHKO.OI.O.  IVW  ■  MIIKM1.?) 

3?  I  b  I  AI  AM  |  KK  )  .  *-  "  .MS  I  M  )]  Ml  *0 
GU  TU  30 
*0  K  •  *♦ 1 

ATAIK  I  •  AT 

tU  I K  )  *  Ml  -  KIM, 

1 1- 1  A  A  .  t  J.CiiDk  I  I )  Iks  I  Jl  »tt  I  K  J 
GO  TU  30 

C  t  A  I  AMt  I  Am  H  t  ACM  ItNC.lMS  ASM  NDMMAl.  m  A  I  tM  SltH'tS 
<*1  KMI  J)  »  KK  J  V 

IM  1 1  IN.  to.  1  IKM(J|*KK  <.)»«•>?  MO. 

!*-(  Ids  .to.?  )  kki  j )  »xm  i  j )  •  j  non. 

K  AN  j  AMS  I  <  HOA  T  <  I  Ate  (  J  I  I  /  lOOO.-MI.DA  I  I  I  AfcO  I  J  I  /  1  (100  I  )  •  !  0(10.  I 
1 1 1  J.GI  .1  I  5C<  J  >»  IKSIJ-1  >-MAI.I)  1/  I  KK  ( Jl  -Km  I  J-I I  I 
C  CMfcCK  )-UM  iNOKtAAlNI.  (Wl)tK  DM  A  I  A  I  I  I  IN  A 
013  h  t  «?  .K 

6  IKIAIAI  I  I.LT.AIAI  l-I  ))  MK  I  NT  10m.  Cl  Hit  I  tv  I  ,  CDOt  I  ?0  I  .CIIDt  IIM 
1  t  I  I  MM  .60. OM. I)  TO  <•? 

XKM  a  KMT 
DDT  l»?.KK 

?  t H A  I  AM  I  |  )  . L  I  . A  I  AM  I  I  -  1  I  IKK  IN t  1 OM ,  COUt < ?1  > .CODt I ?? ) .CnOfrl  IN  I 

C  CUNAIKUC.1  *1  CAMUS 
A?  it!  J.r.J.l  Ill  .J-l 

I  *- 1  I  StC  I  .11  .o!  .OK.II  HI  *<. 

ISCN  «  -IStClJ) 

DU  3  I « | • J 

3  I-  |  I  Af.N.MO.  |  AtC  (  I  I  III  •! 

KKP  a  X  X M  I  )1  > 

A*  H- {  J.l.f  .  1  »  At  I.J  I  *1  MAI  J  l-MAI  I  II  I /I  KKI  .1  l-KKI  I  1  I) 
lH  a  0. 

I  •  0 

3«l  •  I  •  | 

XKK I J  |  a  KK ( J ) »0X IK  I  j  I 
it  |  t  .1.1  .1  »  G4I  III  3  V 
IM.I.eo.l  IGU  III  A3 

S'*  HI.I3I  a  KKK(.I)  -  IKK 
*3  Ml.  I  I  I  a  HU  3  I  ♦!)*(.  I  I  I 
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Test  Problem 


4.  As  a  demonstration  of  generating  a  HEC-2  input  deck  from  stream 
cross-section  data,  a  creek  sector  wa3  hypothesized  containing  several 
regular  sections  as  well  as  auxiliary  or  transfer  sections,  three 
bridges,  a  tributary,  and  several  locations  with  changing  n  values.  The 
creek  sector  was  assumed  to  be  located  in  minibasins  12  to  14  of  sub¬ 
basin  •;(),  as  shown  in  Figure  A3. 

t.  Crc <s-section  40121  is  a  regular  ground  section  with  two 
auxiliary  sections  located  1000  and  500  ft  downstream  from  the  main 
section  and  one  auxiliary  section  4 00  ft  upstream.  The  downstream  sec¬ 
tion  start i  with  three  r  ighness  values  on  an  NC  card  superceded  by  a 
set  of  five  new  roughness  values  at  the  main  section  to  be  punched  on 
an  fiH  card. 

6.  Section  40124  cuts  along  the  center  line  of  a  bridge  designated 
for  the  Normal  Briige  Routine.  Downstream  and  upstream  auxiliary  sec¬ 
tions  are  located  a*  distances  -200  and  100  ft  from  the  main  section, 
respectively.  A  new  NC  card  should  be  entered  at  the  main  section,  then 
changed  back  to  the  old  NC  card  upstream  fr  m  the  bridge.  X3  cards  for 
artificial  levees  are  called  for  i  vnstream  and  upstream  from  the  bridge. 

r.  Section  40128  is  a  regular  ground  section  200  ft  downstream 
from  the  inaccessible  Junction  with  a  tributary.  An  auxiliary  section 
is  to  be  placed  20  ft  above  the  Junction  and  labeled  40129. 

Section  40132  is  a  "Special  Bridge,"  located  1020  ft  upstream 
from  the  Junction.  An  auxiliary  section  is  to  be  placed  1000  ft  down¬ 
stream  from  the  main  section  and  two  more  auxiliary  sections  are  needed 
100  ft  upstream  and  downstream  from  the  main  section.  Again,  changes  in 
NC  cards  at  'he  bridge  and  X3  cards  will  be  needed.  The  bridge  has  a 
weir  flow  over  a  roadbed  of  variable  profile;  therefore,  the  limits  and 
profile  of  the  weir  are  to  be  entered  onto  BT  cards. 

9.  Section  -40129  is  a  section  on  the  tributary  immediately  above 
the  Juncti  n.  The  surveyed  section  is  located  70  ft  upstream  from  the 
Junction,  and  an  auxiliary  section  is  to  be  established  50  ft  downstream 
from  the  surveyed  section  and  is  to  be  labeled  -40129. 
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Figure  A3.  Test  creek  for  XSEC-PROC  demonstration 


10.  Section  U01U2  is  a  special  bridge  section  with  weir  flow  over 
a  horizontal  profile  length  of  300  ft.  This  length  will  have  to  be 
indicated  on  field  t  of  the  SB  card,  and  the  low  chord  as  well  as  top  of 
the  road,  including  railing,  are  to  be  entered  into  fields  h  and  5  of 
the  X2  card  following  the  first  XI  cards  above  the  bridge.  Two  auxil¬ 
iary  sections  are  to  be  placed  100  ft  downstream  and  upstream  from  the 
bridge. 

11.  Tables  A1  and  A 2  are  listings  of  the  XSEC-PR0C  input  and  out¬ 
put  deck3.  It  can  be  seen  that  the  output  constitutes  a  complete  HEC-2 
input  deck  with  the  exception  of  T  and  J  cards. 
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0.0 

0.0 

0.0 

0.0 

0.0 

718.0 

717.) 

n 

147.4 

76. 

70. 

70. 

1.00 

0.00 

t?o 

0160 

770. 


9S0. 

759. 

700. 

0.0 

wo. 

0.0*0 

wo. 

boo. 

165. 

JOO. 

0.0 

714.0 


?J0. 

m. 

0.0 


APPENDIX  8:  COMPUTER  PROGRAM  HYDPAR2  MANUAL 


Background 


1.  Computer  program  KYDPAR2  ia  an  amplification  of  the  HYDPAR  com¬ 
puter  program  developed  by  the  Hydrologic  Engineering  Center  (HEC). 

Vher eas  HYDPAR  ia  essentially  a  bookkeeping  program  which  computes  many 
of  the  input  values  for  HBC-1,  HYDPAF2  Is  designed  to  produce,  in  card 
form,  90  to  95  percent  of  the  input  to  HBC-1,  in  v.  e  specified  sequence, 

as  well  as  QT  and  J1*  cards  for  the  preliminary  run  of  HEC-2. 

# 

2.  HTDPAR2  is  written  to  distinguish  between  subbasins  and  mini¬ 
basins  and  identifies  the  outflow  section  of  a  minlbasin  by  a  ‘3-digit 
number.  The  first  two  digits  identify  the  subbasin,  the  third  and  fourth 
digit  the  minlbasin,  and  the  fifth  digit  is  reserved  for  cross  sections 
within  the  minibasins.  Subbasins  and  minibasins  should  be  numbered  in 
upstream  ( HBC-J !  sequence, 
and  the  section  number  at 
the  minlbasin  outlet  is 
zero.  Minlbasin  inlet 
points,  with  section  num¬ 
ber  9,  are  identified  by 
HYDPAR2  only  when  this 
inlet  is  a  Junction  of 
flow  from  two  or  more 
ainiba3ins.  For  example, 
if  the  subbasin  shown  in 
Figure  81  is  labeled  No. 

19,  HYDPAR2  will  compute 
drainage  areas  at  tne  in¬ 
let  and  outlet  points  of 
minlbasin  2,  labeling  the 
locations  19029  and  19020, 
respectively.  Since  minl¬ 
basin  1  receives  inflow 


Figure  B1 .  Example  watershed 
with  12  aubbasins 
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only  from  minibasin  2,  no  section  19019  is  required. 

3.  The  program  can  handle  up  to  5  subbasins  with  up  to  50  mini¬ 
basins  each.  It  will  scan  the  master  data  file  for  grid  areas  belonging 
to  any  one  of  the  designated  subbasins,  store  all  pertinent  data  in  its 
memory,  then  proceed  to  process  the  data  from  each  subbasin  separately 
into  printed  and/or  punched  output,  as  described  below. 

U.  H YD PAR 2  can  perform  two  routines,  labeled  1RD  =  1  and  2,  respec¬ 
tively.  In  IRD  *  1 ,  the  program  adds  up  grid  points  belonging  to  all 

designated  minibasins  and  determines  not  only  the  areas  of  the  ml  r.i- 
b as ins ,  but  also  the  drainage  area  contributing  to  each  minibasin  outlet. 
In  Figure  81,  for  example,  the  drainage  basin  computed  for  minibasin  3 
includes  the  areas  Tor  minibasins  3,  * ,  5»  6,  and  7.  From  these  drainage 
areas  a  set  of  preliminary  discharges  are  computed  by  the  formula 

q,  -  CJi#5(DA)0-7  (1) 

where 

Q  »  preliminary  discharges 

J 

C  »  an  arbitrary  constant  to  be  chosen  on  the  basis  of  the  desired 
range  of  flows 

J  *  1,  2....NSQ,  a  counter  which  produces  N5Q  flow  rates  up  to  a 
maximum  of  nine  flows 

DA  »  disposal  urea,  acres 

The  1.5  exponent  creates  a  parabolically  increasing  set  of  flow  rates, 
allowing  a  high  upper  bound.  In  the  Wolf  River  Project,  C  was  made 
equal  to  2.0.  With  these  values,  for  example,  a  2000-acre  drainage 
basin  can  be  assigned  preliminary  discharges  varying  between  UQO  and 
10.&00  cfs  roughly  corresponding  to  0.2  and  5  in. /hr  of  runoff.  These 
preliminary  flows  are  used  by  HEC-2  to  produce  a  storage-outflow  table 
for  each  minibasin.  In  addition  to  the  QT  cards  which  dictate  the  dis¬ 
charges  used  by  RBC-2,  HYDPAR  also  produces  the  J U  cards  which  designate 
all  minibasins  for  which  storage  outflow  tables  are  requested. 

5.  In  routine  IRD  a  2,  the  Input  variables  for  HB2-1  are  read  from 
instruction  cards  or  computed  and  punched  out  in  the  sequence  required 
by  HEC-1 .  As  one  possible  option,  rainfall  cards  comprising  a  design 


hyetograph  can  be  read  and  multiplied  by  up  to  10  factors  to  produce 
multiples  of  this  hyetograph.  From  land  use,  hydrologic  soil  group  and 
land  slopes  read  in  the  master  data  file,  weighted  average  Soil  Conserva¬ 
tion  Service  (SCS)  curve  numbers,  and  unit  hydrograph  lag  times  are 
computed. 

6.  The  hydrologic  effect  of  CN,  namely  runoff  or  runoff  peak,  is 
not  proportional  to  CN.  Although  commonly  done,  the  arithmetic  average 
curve  number,  CN ,  can  be  shown  to  underestimate  runoff,  especially  for 
small  rainfalls.  A  more  acceptable  and  accurate  method  of  averaging 
the  CN  values  of  grids  within  a  runoff  area  {minibasin,  in  this  case) 
should  involve  weighting  on  the  basis  of  runoff  effects.  A  routine  to 
accomplish  this  was  incorporated  in  HYDPAR2. 

7.  According  to  the  SCS  formulas: 


S 


1000 

CN 


10 


(2) 


and 

2 

(If  P«xS,  Q*0 ) 

which,  through  an  exchange  of  terms,  becomes 


S 


*S) 


(3) 


where 

S  ■  storage  capacity  of  the  soil  (inches). 

P  ■  cumulative  rainfall  since  start  of  event  (inches). 

Q  ”  cumulative  runoff  since  start  of  event  (inches). 

x  •  portion  of  storage  capacity  S  to  be  used  as  initial  abstrac¬ 
tion  of  rainfall  P  (SCS  uses  x  *  0.2). 

8.  Procedures  for  computing  runoff-weighted  CK  values  are  as 
follows.  For  a  runoff  area  containing  n  values  of  CN  from  n  grids  of 
equal  area  and  the  approximate  total  rainfall  P  ,  4  is  determined 
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for  each  of  the  n  grids  and  an  arithmetic  average  value,  Q  ,  computed. 

A  runoff- weighted  curve  number,  CN*,  is  then  determined  through  an 
iterative  solution  of  Equation  3.  The  iterative  process  is  initiated 
by  substituting  ■  1000/CN  -  10  on  the  right  side  and  solving  for 

on  the  left  side.  If  |S0  -  /S^  j  <  0.01  ,  5^  iB  accepted  as  the 

solution.  If  not,  a  new  is  set  equal  to  S0  ♦  S^/2  and  solution 

of  Equation  3  is  repeated.  When  an  acceptable  S  has  been  found,  a 
runoff-weighted  average  curve  number  is  computed,  i.e.,  CM*  =  1000/S+10. 

9.  A  routine  was  added  to  recode  the  SCS  erosion-slope  index, 
listed  in  column  6  of  the  master  data  file  into  separate  erosion  and 
slope  parameters  for  HEC  input.  The  SCS  erosion-slope  index  is  con¬ 
verted  to  individual  indices  by  the  following  code: 

a.  The  first  digit  is  the  erosion  index. 

b.  The  last  digit  is  the  slope  index,  except  when  the  erosion 
index  is  greater  than  2,  in  which  cane  the  slop*  index  is 
the  last  digit  plus  1. 

c.  The  slope  index  is  further  converted  to  a  percent  slope 
according  to  the  following  code: 

Slope  index:  1  2  3  fc  5  6  7 

Percent  Slope:  1.0  U.O  8.0  12.5  20.0  35.0  8.0 

Example : 

SCS  Erosion-  Erosion  Slope 

Slope  Index  Index  Index  Slope  Percent 

201  21  1 
302  3  3  8 

d.  The  input  to  the  HEC-1  K  cards  is  produced,  which  provides 
the  operational  direction  in  the  HEC-1  program.  To  obtain 
these  K-cards  in  the  proper  sequence,  the  following  rules 
should  be  observed: 

(1)  Starting  with  the  longest  streambed,  number  the  sub¬ 
basins  consecutively  in  upstream  direction  (l  to  5  in 
Figure  Bl). 

(2)  Pick  up  the  most  upstream  tributary  and  label  the 
basins  again  in  upstream  order  (6  to  7  in  Figure  Bl). 

(3)  Repeat  step  2  with  the  next  lower  tributary,  numbering 
the  basins  as  before  (8  to  10  in  Figure  Bl).  If  the 
tributary  has  subtributaries  itself,  continue 


numbering  these  basins  on  the  subtributaries  (11  and 
12  in  Figure  Bl)  before  moving  down  to  the  next 
tributary  of  the  main  stream. 

Input  Directions 


10.  For  run  1,  the  input  must  consist  of  one  A1 ,  J1 ,  J2,  and  J3 
card  each,  followed  by  a n  A2,  H,  and  one  or  more  K  cards  for  each 
subbasin. 

11.  For  run  2,  the  full  sequence  of  input  cards  shown  on  the  fol¬ 
lowing  pages  (which  include  all  of  the  run  1  input  cards)  must  be 
entered. 

1. .  Values  may  be  entered  in  I  or  F  type  except  where  I  type  is 
specifically  indicated  under  "Value." 

AX  Card  Title  Card,  Format  (2X,  Ait,  11A6) 

1  i.  A1  may  be  punched  (optionally)  in  columns  1,  2,  followed  by  a 
title  of  alphanumeric  type  in  columns  3-72. 

J1  card  Operational  Instructions,  Format  (2X,  FG.O,  8F8.0) 


Field 

Columns 

Variable 

V  alue 

0 

1-2 

J1 

1 

3-8 

KSUB 

♦ 

2 

9-1 6 

NOISE 

♦ 

3 

17-2U 

Nsq 

♦ 

U 

25-32 

SIZE 

9.87 

5 

33-<*0 

ILURFT 

- 

0  or 

6 

bl  — 1*8 

IP 

0 

1 

7 

1*9-56 

IRD 

1 

2 

8 

57-6U 

LLN 

0 

_ Description _ 

Card  Code 

Number  of  subbasins  (max  5) 

Number  of  land  uses  (max  25) 

Number  of  storage  -Q  values  (max  9) 

Grid  cell  sire  in  acres 

No  printout  of  minibasin  sizes, 
slopes,  CN  values,  etc.,  is 
wanted 

Normal  printout  requested 

No  punched  cards  requested 

launched  cards  from  run  1  ( RFC -2  QT 
cards)  or  from  run  2  ( HEC-1  in¬ 
put)  requested 

Run  1  requested 

Run  2  requested 

CN  values  will  be  computed  by 
H  YD  PAR 
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J1  'ard  (Continued) 


Field 

Columns 

Variable 

Value 

Description 

1 

CN  values  are  read  from  Data  Base 

9 

65-72 

NRNG 

0 

No  rainfall  will  be  read 

♦ 

No.  of  Hyetographs  to  be  read 

10 

73-80 

NXRN 

0 

If  NRNG  -  0 

♦ 

No.  of  multipliers  to  be  read  in 

R1  card.  Also  No.  of  new 
hyetographs  to  be  created 

J2  Card 

Data  Base  Access, 

Format  (2X 

,  F6.0,  8F8.0) 

lk 

.  The  values  in  the  J 2  card  are  the  fields  in  the  Data  Base 

under  which  the 

indicated 

variables  can  be  found. 

Field 

Columns 

Variable 

Value 

Description 

0 

1-2 

J2 

Card  Code 

1 

3-8 

ISD 

♦ 

Subbasin  number 

2 

9-16 

IBS 

♦ 

Minibasin  number 

3 

17-2k 

ISL 

♦ 

Soil  group 

k 

25-32 

IHYO 

♦ 

Hydrologic  soil  group 

5 

33-kO 

1SLP 

♦ 

Slope-erosion  index 

6 

kl-U8 

ILD 

♦ 

Land  use  category 

7 

1*9-56 

ICS 

0 

No  CN  values  will  be  read  from 
data  base 

♦ 

Data  base  field  where  CN  values 
are  listed 

J3  Card 

Format 

(2X,  F6.0, 

F8.9,  518) 

Field 

Columns 

Variable 

Value 

Description 

0 

1-2 

J3 

Card  Code 

1 

3-8 

XQ 

♦ 

Coefficient  for  preliminary  design 
flows  (see  Equation  1) 

2 

9-16 

P 

♦ 

Precipitation  in  inches  used  in 
weighted  CN  averaging;  recom¬ 
mended  value  *  3 

3 

17-2k 

ISUB(l) 

I 

First  subbasin  number 

k 

25-32 

ISUB(2) 

I 

Second  subbaaln  number 

J3  Card  (Continued) 


Field 

Columns 

Variable 

Value 

Description 

5-7 

Repeat  for  total 

of  NSUB,  given  in  J1  card 

8 

57-6U 

1RTIMF 

I 

Rainfall  and  routing  time  interval 
in  minutes 

9 

65-72 

NRNI 

No.  of  rainfall  intervals  in  each 
hyetograph 

R1  Card 

Only  used  when  IRD 

c  2  ( run 

2)  and  NRNG  >  0  (Jl-9) 

Field 

Columns 

Variable 

Value 

Description 

0 

1-2 

R1 

Code 

1 

3-8 

XRN(l) 

F 

First  Multiplier  to  create  new 
hyetograph 

2 

9-16 

XRN( 2) 

Second  multiplier 

etc . 

repeat 

NXHN  times 

(Jl-10) 

RC  Cards  Only  used  when  R1 

cards  are  used 

Field 

Columns 

Variable 

Value 

Description 

0 

1 

2 

Card  Code.  (Note  that  only  col- 

umn  1  is  used  for  the  code,  in 
order  to  allow  the  use  of  the 
sam>-  rainfall  data  cards  as 
input  to  program  KEC-l) 


1 

2-8 

PR EC  ( 1 ) 

F 

Rainfall,  in  units  to  be 
HEC-1  (usually  inches) 
during  first  interval. 

used  in 
failing 

2 

7-16 

PREC(2) 

F 

Rainfall  in  second  interval 

etc. , 

repeat 

NRNI  times 

( J3-9) 

DU  Cards 

Land 

Uses,  Format 

(1*A6,  LF8.0) 

15. 

Only  used  when  IRD 

a  2  (run 

1)  Repeat  NLUSE  times  as 

given 

on  J1  card. 

Field 

Columns 

Variable 

Value 

Description 

1-3 

1-2U 

LU 

A-N 

Land  use  name 

1* 

25-32 

CVNM(l) 

♦ 

CN  for  hyd.  group  A 

5 

33-1*0 

CVNM(2) 

♦ 

CN  for  hyd.  group  B 

6 

1*1— 1*8 

CVNM(3) 

♦ 

CN  for  hyd.  group  C 

7 

1*9-56 

CVNM(1») 

♦ 

CN  for  hyd.  group  D 
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16.  The  following  set  of  cards,  A 2  to  K,  will  have  to  be  repeated 
for  each  subbasin. 

A2  Card  Title  Card  for  subbasin  described  same  format  as  kl  card 


HI  Card 

Specific 

Instruction 

for  HEC- 

-1  Input,  Format  (2X,  12,  9I1*, 

UF4.0) 

Field 

Columns 

Variable 

Value 

Description 

0 

1-2 

HI 

Card  Code 

1 

3-4 

NMIti 

I 

Number  of  minibasins 

2 

5-8 

NRR 

I 

Number  of  branches 

3 

9-12 

NJ4 

I 

Number  of  sections  for  which 
storage-outflow  cards  are 
generated  by  HEC-2 

4 

13-16 

HI 

I 

Entry  for  HBC-1  card  M,  field  1 

5 

17-20 

M2 

I 

Entry  for  HEC-1  card  M,  field  2 

M2  «  2  if  SCS  hydrographs  are 
generated 

6 

21-24 

I<U 

I 

Entry  for  HBC-1  card  Q,  field  1 

0 

If  HI  i  2 

Number  of  nonrecording  rain  gages 
if  weighting  factors  are  con¬ 
stant  over  the  subbasin 

♦ 

Number  of  nonrecording  rain  gages 
if  weighting  factors  must  be 
read  for  all  minibasinB 

7 

25-28 

102 

I 

Entry  for  HEC-1  card  Q,  field  2 
Same  as  Q1 ,  but  for  recording 
rain  gages 

3 

29-32 

JY4 

I 

Entry  for  HEC-1  card  Y,  field  4 

9 

33-36 

JY11 

I 

Entry  for  HEC-1  card  Yl,  field  1 

10 

37-40 

JY17 

I 

Entry  for  HBC-1  card  Yl,  field  7 

11 

41-44 

X2 

F 

Entry  for  HEC-1  card  X,  field  2 

12 

45-48 

X3 

F 

Entry  for  HBC-1  card  X,  field  3 

13 

49-52 

Y1 

F 

Entry  for  HBC-1  card  Y,  field  1 

14 

53-56 

Y2 

F 

Entry  for  HBC-1  card  Y,  field  2 

J‘»  Cards  Locations  for  Storage-Outflow  Data 


IT.  This  is  the  same  card,  or  cards,  which  are  submitted  as  HEC-2 
input.  The  first  three  fields  on  the  first  JU  card  contain  data  spec¬ 
ified  by  HEC-2,  and  the  following  fields  in  a  10K8.0  format  contain 
NJU  section  numbers,  on  as  many  cards  as  needed. 

Storage-Outflow  Cards 

18.  When  NJU  section  numbers  are  specified  on  the  Jl*  cards,  HEC-2 
will  produce  storages  for  all  reaches  between  successive  sections  and 
outflows  at  the  corresponding  downstream  sections.  Together  with  the 
Y  and  Y1  cards,  which  are  not  related  to  the  storages  and  outflows  on 
the  Y2  and  Y3  cards,  HEC-2  will  thus  cause  (NJl-l)  set3  of  four  cards  to 
be  punched  out.  This  entire  deck  is  entered  into  HYDPAR2  after  the  Jl* 
cards.  HYDPAB2  will  then  reproduce  the  storage-output  cards,  but  enter 
zero  values  in  field  1,  and  place  the  cards  into  the  proper  place  in 
the  HEC-1  input  deck. 

J 5  Card ( s )  Only  use  if  I»il  <  0,  indicating  nonrecording  gage3  with  con¬ 
stant  weighting  factors.  Format  (.'X,  12,  Flt.O,  7(11*,  Fl.O)) 


Field  Columns  Variable  Value 


JDNR(l) 
VTSH  ( 1 ) 
JDNR(2) 
V7SB ( 2 ) 


Description 


Card  Code 


13-16 


I.D.  So.  of  first  nonrec.  gage 
Weighting  factor  for  gage 

Same  as  above,  for  second  gage 


Continue  for  1 01  g&geB. 


J 6  Card  Only  used  if  IQ2  <  0 

10.  Sam#*  as  J5  card,  but  for  recording  gages. 
J7  Card  Only  if  I(U  >  0  (Variable  weighting) 


Field  Columns  Variable  Value 


IDNP(l) 


Description 


Card  Code 


I.D.  number  of  nonrec.  rain  gage 


blank 
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Unlabeled  cards  following  J7  card 


Field  Columns  Variable  Value 

1  1-1*  WTNR(1,I)  F 


Description 


Weighting  factor  for  gage  IDNR(l) 
applicable  on  minibasin  1 

Weighting  factor  for  some  gage 
but  minibasin  2 


etc,,  for  all  minibasina,  even  if  weighting  factor  is  zero. 

Repeat  J 7  and  weighting  factor  cards  for  all  nonrecording  gages. 

Cards  J8  Only  if  IQ2  >  0  (Variable  weighting) 

20.  Same  as  J7  card(s),  but  for  recording  gages. 

L  Cards  Flow  Path  Lengths,  Format  (2X,  P6.0,  9F8.0) 

21.  These  cards  provide  the  flow  path  lengths  for  all  minibasins 
within  the  subbasin  under  consideration.  These  lengths  are  used  to 
compute  the  lag  time  for  the  unit  hydrograph  derivation  in  each  mini¬ 
basin  and  should  thus  be  estimates  of  the  longest  path  from  the  mini¬ 
basin  outlet  to  any  point  on  the  basin  except  that  narrow  appendices  of 
a  minibasin  should  be  ignored  (see  sketch). 

PC.  For  minibasins  1  to  KM1N,  the 
flow  path  lengths,  in  feet,  should  be  .g 

entered  in  fields  1  to  10  (8  columns  /  J 

each)  on  as  many  cards  as  needed.  7  J  t 

K  Cards  Routing  Instructions  to  pro-  I  *//  / 

duce  HEC-1  K  cards.  Format  (201b)  J  S*//  [ 

23.  These  cards  contain  mini-  /  7 /  l 

basin  numbers  arranged  in  branch  order  S Cf 

required  to  generate  the  proper  se- 
quence  of  K  cards  for  hydrograph,  gen¬ 
eration,  routing,  and  combining  in  HEC-1. 

21*.  Consider  the  stream  with  tributary  in  the  sketch  below.  HEC-1 
will  proceed  generating  and  routing  hydrographs  from  minibasin  b  to  the 
top  of  basin  2,  then  from  basin  5  to  the  top  of  baBin  2,  combine  hydro¬ 
graphs  and  continue  down  to  the  bottom  of  basin  1  (or  top  of  basin  0). 


APPENDIX  C:  COMPUTER  PROGRAM  RAINWEIGHT 


Weighting  Techniques 


1.  Typical  conventional  rainfall  weighting  techniques  are  averag¬ 
ing,  Thiessen,  and  Isohyet&l  contours. 

2.  Recently,  it  has  become  popular  to  weight  rainfall  readings  at 


neighboring  stations  by  the  inverse  of  the  distance  to  any  point  under 
consideration.  If  we  have,  for  example,  three  rainfall  readings  R  at 


gages  1  to  3  and  want  a  weighted  rainfall  at  point  P  this  rainfall  is 


where  is  the  distance  from  point  P  to  gage  i.  If  there  are  more 

than  three  gages,  often  only  the  three  closest  gages  are  considered 
( Pi gure  Cl ) . 

3.  If  a  weighted  average  rainfall  is  wanted  for  an  entire  water 
shed,  the  watershed  is  divided  into  a  grid  network  for  at  least  50 
squares,  as  shown  in  Figure  C2. 


Figure  Cl.  Watershed  with 
three  rain  gages 


Figure  C2.  Watershed  divided 
in  grid  network 
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h.  EhcH  grid  area,  including  the  boun  iary  areas  larger  than  one 
half  square,  is  given  a  set  of  coordinates,  (fc,i>)  for  example  Tor  the 
cross-natched  area.  The  rain  gages  are  given  the  coordinates  of  the 
square  they  are  located  in.  Now  the  distance  from  any  grid  area  to  a 
neighboring  gage  is  given  by  yjux  )2  ♦  (Ay)^  .  The  distance  from  the 
shaded  area  to  gage  1,  for  example,  is  ♦  2^  *  2.83  . 

5.  The  weighted  rainfall  for  each  grid  space  is  then  computed  by 
Equation  1,  and  the  weighted  rainfall  arxjur.t a  for  all  grid  spaces  are 
averaged.  Any  grid  space  in  which  a  gage  is  located  in  automatically 
weighted  100  percent  in  favor  of  that  gage,  and  all  other  gages  are 
ignored. 

3* 

6.  Dean  and  Snyder  suggest  using  the  inverse  of  distance  squared 
(i.e.,  k  ■  2)  as  a  weighting  factor,  but  differences  in  resulting 
weighted  rainfalls  are  usually  small. 

7.  The  procedure  is  obviously  very  tedious  unless  a  computer  pro¬ 
gram  is  used. 


8.  For  weighted  rainfall  computing  by  program  RAINVKIGHT,  it  is 
nly  necessary  to  *•,<•«!, fy  the  number  f  rain  gages,  their  co  rdi nates, 
the  number  of  grid  c  lunns,  and  the  range  of  rows  covered  by  each  grid 
column.  The  input  directions  are  as  follows: 

Card  1:  NO,  NC,  IFR,  XN  Format  (  3I*« .  FU .  1 ) 

Card  2:  X(l),  Y(l),  X(2),  Y(2)  _  X(NG),  Y(NG)  Format  (20F14.0) 

Card  3:  YT(1),  YB(1).  YT(2),  YB(2) . YT(NC).  YB(NC) 

Format  (20FU.0) 

Card  !*:  PR(l),  PR(2)  _  PR  (KG)  Format  (20F4.1) 


NG  »  Number  of  gages 

NC  »  Number  of  grid  columns 


See  References  at  the  end  of  the  main  text 


IPR  ■  Printing  instruction.  Set  IPR  *  1  if  »  printout 
is  wanted  of  the  three  gages  nearest  to  any  grid 
space. 

XN  "  Exponent  of  distance  to  be  used  in  weighting. 
X(i),  Y(i)  =  x  and  y  coordinates  of  gage  i 
YT( J  ) ,  YB(j)  *  'i’op  and  bottom  row  number  in  column  1 
PR( i )  ■  Rainfall  reading  at  gage  i 
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APFKNDIX  D:  REFERENCE  MATERIAL  ON  OPTIMAL 
ENERGY  SLOPE  AVERAGING  AND  REACH  LENGTH 


1.  One  of  the  major  sources  of  error  in  open  channel  backwater  com¬ 
putations  lie3  in  the  determination  of  average  energy  gradients  over 
relatively  Ion*?  reaches  of  river  channel.  The  two  questions  the  hydrau¬ 
lic  engineer  is  confronted  with  in  this  regard  are: 

a.  What  reach  lengths  should  be  used? 

b .  Given  the  down at  rear,  and  upstream  energy  slopes  within  a 
reach,  how  should  these  slopes  be  averaged  to  determine 
head  losses? 


2,  In  all  likelihood,  the  accuracy  of  backwater  computations  in¬ 
creases  with  increasing  number  of  cross  sections  or  decreasing  reach 
length.  The  cost-conscious  engineer  must  keep  in  mind,  however,  that 
the  surveying  of  cross  sections  is  a  fairly  costly  task  and  that  it  is, 
therefore,  desirable  to  moxir.iie  the  reach  lengths  in  a  channel  without 
Introducing  excessive  errors  in  the  water  surface  profile  computations. 

i.  Secondly,  most  engineers  have  at  one  time  or  another  seen  sev¬ 
eral  formulas  for  averaging  slopes.  Most  manuals  propose  simple  arith¬ 
metic  averaging  of  upstream  and  downstream  energy  gradients,  and 

Sj  ,  respectively,  according  to  the  formula 


S 


S  ♦ 

u 


2 


(1) 


where  S  is  the  average  energy  gradient  slope.  Others  advocate 
geometric  averaging  by  the  formula 


§  •  Vs u  "  Sd  (2) 

2/3 

or  the  averaging  of  conveyances,  l.l*9/n  AF.'  which  is  equivalent  to 
averaging  the  slopes  taken  to  the  (-1/2)  power  (A  «  cross-sectional 
area;  R  *  hydraulic  radius). 


D1 


It* 

u.  Tavener  presents  some  analytical  answers  to  these  questions. 

In  the  first  place  he  suggests  distinguishing  between  M-l  and  M-2  sur- 
face  profiles,  depending  on  whether  the  energy  slope  is  flatter  or 
steeper,  respectively,  than  the  channel  slope.  He  then  suites  the  inter¬ 
esting  observation  that  the  M-l  curve  resembles  a  parabola  with  a  verti¬ 
cal  axis  while  the  M-2  curve  resembles  a  parabola  with  a  horizontal 


axis.  On  the  basiB  of  this  observation,  Tavener  proceeds  to  show  that, 
for  an  .M-l  curve  the  arithmetic  slope  averaging  equation  (Equation  1)  is 
the  moat  correct  or.e,  while  for  an  M-2  curve  the  formula 


should  be  used. 

5.  Through  a  semi analyt leal  and  experimental  process  Tavener  then 
derived  a  formula  for  the  maximum  allowable  reach  length  L  ,  i.amely 


where 

V  *  Froude  Number,  dimensionless 

S  •  hydraulic  gradient,  ft/ft 

T  *  top  width  of  the  wetted  area,  ft 

2 

A  *  cross-sectional  wetted  area,  ft 

P  ■  wetted  perimeter,  ft 

B  ”  change  in  P  with  reppect  to  depth,  ft /ft 
and  the  subscript  d  denotes  that  the  above  dimensions  correspond  to 
the  downstream  section  of  a  reach. 

6.  By  compari son  of  the  results  from  Equation  U  with  his  experi¬ 
mental  results,  Tavener  found  that  the  equation  slightly  overestimated 
the  allowable  reach  length  ALOVL  and  suggested  the  design  formula 


ALOWL 


(5) 


L 


in  which  c  equals  2.0  for  M-l  curves  and  1.2  for  M-2  curves. 

7.  Equation  5  cannot  be  applied  to  channels  prior  to  the  backwater 
computations  because  all  of  the  dimensions  used  in  the  formula  depend  on 
downstream  depth.  The  equation  is  rather  used  as  a  check  after  each 
backwater  step  to  see  whether  additional  measured  or  interpolated  cross 
sections  should  be  inserted  in  the  following  reach. 

8.  It  may  be  noted  that  Equation  U  is  somewhat  cumbersome  to  per¬ 
form  after  each  backwater  step.  While  A.  ,  P .  ,  and  5 .  are  usually 

u  a  a 

listed  in  the  backwater  tabulations  and  the  top  width  can  be  easily  ob- 

2 

tained,  >\;  and  B,  would  have  to  be  computed  separately.  Further¬ 
more,  ,  the  change  in  P  with  depth,  could  assume  unreasonably 

high  values  if  the  bank  slope  changes  abruptly  with  depth.  By  testing 
a  large  range  of  parabolic  cross  sections,  the  value  of  2B/3P  was 
found  to  vary  between  12  and  20  percent  of  5T/3A  .  The  writer,  there¬ 
fore,  decided  to  climate  the  term  2/3  B/P  by  setting  3/2  B/P  =  0.2 
5/3  T/A  and  to  replace  r  by  Q‘T/gA  .  With  these  substitutions 
Equation  5  is  changed  to 


2 

in  which  g  is  the  gravitational  constant  (32.2  ft  /sec*' )  and  c  equals 
1.5  for  M-l  curves  and  1.0  for  M-2  curves. 

9.  Equation  6  looks  even  more  complicated  than  Equation  *»,  but 
lends  itself  very  well  to  nomograph  solution.  The  nomograph  shown  in 
Figure  1  was  constructed  and  is  used  as  follows: 

a.  Find  Q/A  and  A/T  for  the  downstream  section  of  a  new 
reach  and  draw  line  1  from  scale  2  through  scale  3  to 
scale  ka,  where  the  f2  is  read. 

b.  Subtract  V*  from  unity  and  draw  a  line  2  from  1-F*  on 
scale  bb  to  A/T  on  scale  2. 
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EXAMPLE  I  f  ROM  TAVENER 
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Figure  Dl.  Nomograph  for  determining  maximum 
allowable  reach  length*4 
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n 


for  each  of  the  n  grids  and  an  arithmetic  average  value,  Q  ,  computed. 

A  runof f-weighted  curve  number,  CN*,  is  then  determined  through  an 
iterative  solution  of  Equation  3.  The  iterative  process  is  initiated 
by  substituting  =  (1000/CN)  -  10  on  the  right  side  and  solving  for 

on  the  left  side.  If  | (S^  “  S ^ ) /S ^  |  <  0.01  ,  is  accepted  as  the 

solution.  If  not,  a  new  is  set  equal  to  (S^  +  Sj)/2  and  solution 

of  Equation  3  is  repeated.  When  an  acceptable  S  has  been  found,  a 
runof f-weighted  average  curve  number  is  computed,  i.e.,  CN*  =  1000/ (S+10). 

9.  A  routine  was  added  to  recode  the  SCS  erosion-slope  index, 
listed  in  column  6  of  the  master  data  file,  into  separate  erosion  and 
slope  parameters  for  HEC  input.  The  SCS  erosion-slope  index  is  con¬ 
verted  to  individual  indices  by  the  following  code: 

a.  The  first  digit  is  the  erosion  index. 

b.  The  last  digit  is  the  slope  index,  except  when  the  erosion 
index  is  greater  than  2,  in  which  case  the  slope  index  is 
the  last  digit  plus  1 . 

c.  The  slope  index  is  further  converted  to  a  percent  slope 
according  to  the  following  code: 


Slope  index: 

1  2  3 

4  5 

6  7 

Percent  Slope: 

1.0  4.0  8.0 

12.5  20.0 

35.0  8.0 

Example: 

SCS  Erosion- 

Erosion 

Slope 

Slope  Index 

Index 

Index 

Slope  Percent 

201 

2 

1 

1 

302 

3 

3 

8 

The  input  to  the 

HEC-1  K  cards 

is  produced, 

which  provides 

the  operational  direction  in  the  HEC-1  program.  To  obtain 
these  K-cards  in  the  proper  sequence,  the  following  rules 
should  be  observed: 


(1)  Starting  with  the  longest  streambed,  number  the  sub¬ 
basins  consecutively  in  upstream  direction  (1  to  5  in 
Figure  Bl). 

(2)  Pick  up  the  most  upstream  tributary  and  label  the 
basins  again  in  upstream  order  (6  to  7  in  Figure  Bl). 

(3)  Repeat  step  2  with  the  next  lower  tributary,  numbering 
the  basins  as  before  (8  to  10  in  Figure  Bl).  If  the 
tributary  has  subtributaries  itself,  continue 


